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ABSTRACT 


Predicting the liquefaction resistance of soil is an important step in the engineering 
design of new and the retrofit of existing structures in earthquake-prone regions. The procedure 
currently used in the U.S. and throughout much of the world to predict liquefaction resistance is 
termed the simplified procedure. This simplified procedure was originally developed by H. B. 
Seed and I. M. Idriss in the late 1960s using blow count from the Standard Penetration Test. 
Small-strain shear wave velocity measurements provide a promising supplement and in some 
cases, where only geophysical measurements are possible, may be the only alternative to the 
penetration-based approach. This report presents guidelines for evaluating liquefaction 
resistance using shear wave velocity measurements. These guidelines were written in 
cooperation with industry, researchers and practitioners, and evolved from workshops in 1996 
and 1998 as well as review comments received on an earlier draft. The guidelines present a 
recommended procedure, which follows the general format of the penetration-based simplified 
procedure. The proposed procedure has been validated through case history data from more 
than 20 earthquakes and 70 measurement sites in soils ranging from clean fine sand to sandy 
gravel with cobbles to profiles including silty clay layers. Deterministic liquefaction resistance 
curves were established by applying a modified relationship between the shear wave velocity and 
cyclic stress ratio for the constant average cyclic shear strain suggested by R. Dobry. These 
curves correctly predict moderate to high liquefaction potential for over 95 % of the liquefaction 
case histories, and are shown to be consistent with the penetration-based curves in sandy soils. 
From logistic regression and Bayesian models, the recommended deterministic curve is 
characterized with a probability of liquefaction of about 26 %. To further validate the 
procedure, additional case histories are needed with all soil types that have and have not 
liquefied, particularly from deeper deposits (depth > 8 m) and from denser soils (shear wave 
velocity > 200 m/s) shaken by stronger ground motions (peak ground acceleration > 0.4 g). The 
guidelines serve as a resource document for practitioners and researchers involved in evaluating 
soil liquefaction resistance. 


KEYWORDS: building technology; earthquakes; in situ measurements; seismic testing; shear 
wave velocity; soil liquefaction 
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CHAPTER 1 


INTRODUCTION 


1.1 BACKGROUND 


A major cause of damage from earthquakes is liquefaction-induced ground failure. For 
example, direct property loss caused by liquefaction during the 1989 Loma Prieta, California 
earthquake (moment magnitude, M4, = 7.0) was over $100 million (Holzer, 1998). Large 
indirect property loss by fire almost occurred in 1989 when liquefaction-induced ground 
deformation ruptured water mains that served the Marina District of San Francisco. 
Fortunately, the fire in the Marina District at Divisadero and Beach Streets was contained to the 
three-story apartment building where it ignited. It was also fortunate that the 1989 earthquake 
did not occur closer to the San Francisco Bay area. The cities of Kobe and Osaka, Japan were 
not so fortunate. The 1995 Hyogoken-Nanbu earthquake (MM, = 6.9) directly struck this 
metropolitan area, causing over $100 billion in property damage (Kimura, 1996). A significant 
portion of the damage in Kobe can be attributed to liquefaction-induced ground deformation. 
These are just two of many examples of major damage caused by liquefaction-induced ground 
failure. Predicting soil liquefaction resistance is an important step in the engineering design of 
new and the retrofit of existing structures in earthquake-prone regions. 


The procedure widely used in the United States and throughout much of the world for 
predicting the liquefaction resistance of soils is termed the simplified procedure. This simplified 
procedure was originally developed by Seed and Idriss (1971) using blow count from the 
Standard Penetration Test (SPT) correlated with a parameter called cyclic stress ratio that 
represents the seismic loading on the soil. Since 1971, the procedure has been revised and 
updated (Seed, 1979; Seed and Idriss, 1982; Seed et al., 1983; Seed et al., 1985). Correlations 
based on the Cone Penetration Test (CPT) and shear wave velocity measurements have also 
been developed by various investigators. General reviews of the simplified procedure are 
contained in a report by the National Research Council (1985), a workshop report edited by 
Youd and Idriss (1997), and a journal paper by Youd et al. (2001). 


Small-strain shear wave velocity, Vs, measurements provide a promising supplement and 
in some cases, where only geophysical measurements are possible, may be the only alternative to 
the penetration-based approach. The use of Vs as an index of liquefaction resistance is soundly 
based because both Vs and liquefaction resistance are similarly influenced by many of the same 
factors (e.g., void ratio, state of stress, stress history, and geologic age). 


Some advantages of using Vs are (Dobry et al., 1981; Seed et al., 1983; Stokoe et al., 
1988a; Tokimatsu and Uchida, 1990): (1) Measurements are possible in soils that are hard to 
sample, such as gravelly soils where penetration tests may be unreliable. (2) Measurements can 
be performed on small laboratory specimens, allowing direct comparisons between laboratory 


and field behavior. (3) Vs is a basic mechanical property of soil materials, directly related to 
small-strain shear modulus, Gmax, by: 


Gra = p Vs? (1.1) 
where 


p= the mass density of soil. 


(4) Gmax, OF Vs, 1S in turn a required property in analytical procedures for estimating dynamic 
shearing strain in soil in earthquake site response and soil-structure interaction analyses. (5) Vs 
can be measured by the Spectral-Analysis-of-Surface-Waves (SASW) test method at sites where 
borings may not be permitted, such as capped landfills, and sites that extend for great distances 
where rapid evaluation is required, such as lifelines and large building complexes. 


Three concerns when using Vs to evaluate liquefaction resistance are: (1) Measurements 
are made at small strains, whereas pore-water pressure buildup and liquefaction are medium- to 
high-strain phenomena (Jamiolkowski and Lo Presti, 1990; Teachavorasinskun et al., 1994; Roy 
et al, 1996). This concern can be significant for cemented soils, since small-strain 
measurements are highly sensitive to weak interparticle bonding which is eliminated at medium 
and high strains. It also can be significant in silty soils above the water table where negative 
pore water pressures can increase Vs. (2) No samples are obtained for classification of soils and 
identification of non-liquefiable soft clayey soils. According to the so-called Chinese criteria, 
non-liquefiable clayey soils have clay contents (particles smaller than 5 pm) > 15 %, liquid limits 
> 35 %, or moisture contents < 90 % of the liquid limit (Seed and Idriss, 1982). Andrews and 
Martin (2000) refined this criteria to soils with clay contents (particles smaller than 2 pm) > 10 
% and liquid limits > 32 % (by Casagrade-type percussion apparatus) for non-liquefiable clayey 
soils. (3) Thin, low Vs strata may not be detected if the measurement interval is too large 
(USBR, 1989; Boulanger et al., 1997). . 


In general, borings should always be a part of the field investigation. Surface 
geophysical measurements and cone soundings are often conducted first to help select the best 
locations for borehole sampling and testing. Surface geophysical tests usually involve making 
measurements at several different locations, and provide general, or average, stratigraphy for 
sediments beneath the area tested. The ability of surface geophysical methods to resolve a layer 
at depth depends on the thickness, depth, and continuity of that layer, as well as the test and 
interpretation procedures employed. Cone soundings provide detailed stratigraphy at each test 
location for sediments that can be penetrated. The preferred practice when using Vs 
measurements to evaluate liquefaction resistance is to drill sufficient boreholes and conduct 
sufficient tests to detect and delineate thin liquefiable strata, to identify non-liquefiable clay-rich 
soils, to identify silty soils above the ground water table that might have lower values of Vs 
should the water table rise, and to detect liquefiable weakly cemented soils. 


1.2 PURPOSE 


This report presents guidelines for evaluating liquefaction resistance through shear wave 
velocity measurements. The guidelines are based on an earlier report entitled “Draft Guidelines 
for Evaluating Liquefaction Resistance Using Shear Wave Velocity Measurements and 
Simplified Procedures” by Andnmus et al. (1999), which evolved from two workshops. The first 
workshop was held on January 4-5, 1996 in Salt Lake City, and was sponsored by the National 
Center for Earthquake Engineering Research (NCEER). The second workshop was held on 
August 14-15, 1998 also in Salt Lake City, and was sponsored by the Multidisciplinary Center 
for Earthquake Engineering Research (MCEER, formally NCEER) and the National Science 
Foundation (NSF). These two workshops are herein called the 1996 NCEER Workshop and 
1998 MCEER Workshop. The guidelines present a recommended procedure based on the 
suggestions given at the workshops, as well as review comments received on the draft 
guidelines. The guidelines provide guidance on selecting site variables and correction factors 
that are consistent with the shear-wave-based procedure. 


From the comments received on the earlier draft guidelines, several improvements have 
been made to the guidelines. The major improvements include: (1) Much of the background 
and development information is moved to the appendixes to provide a clearer, more practical 
description of the recommended procedure. (2) A table is added to compare advantages and 
disadvantages of the various in situ Vs test methods for liquefaction assessment. (3) A more 
comprehensive description of the database is given, particularly on the nature of the case 
histories of gravelly soils. (4) The results of a probability study of the Vs-based case history data 
are added to calibrate the recommended liquefaction resistance curve and to compare with the 
results of probability studies of the SPT-based case history data. 


1.3 REPORT OVERVIEW 


Following this introduction, Chapter 2 presents the recommended procedure for 
evaluating liquefaction resistance using Vs. Chapter 3 illustrates the application of the procedure 
using two case studies. And Chapter 4 summarizes the recommended procedure and identifies 
issues that remain to be resolved. 


Eight appendixes are included to assist the reader, and to provide information used in the 
development of the guidelines. Appendix A presents a list of references cited in the guidelines. 
Appendix B provides a list of Symbols and Notation, and Appendix C provides a Glossary of 
Terms. Appendix D reviews six proposed Vs-based liquefaction resistance curves. Appendix E 
describes the general characteristics of case history data used to develop the recommended 
liquefaction resistance curves. Appendix F presents the development of the recommended 
curves. Appendix G considers three probability models for the case history data. Finally, 
Appendix H presents a summary of the case history data. 


CHAPTER 2 


LIQUEFACTION EVALUATION PROCEDURE 


This chapter presents guidelines for using the Vs-based liquefaction evaluation procedure 
originally proposed by Andrus and Stokoe (1997) and subsequently updated in the report by 
Andrus et al. (1999) and the paper by Andrus and Stokoe (2000). The evaluation procedure 
follows the general format of the Seed-Idriss simplified procedure, and the general 
recommendations of the 1996 NCEER Workshop (Youd et al., 1997) and 1998 MCEER 
Workshop (Youd et al., 2001). It requires the calculation of three parameters: (1) the level of 
cyclic loading on the soil caused by the earthquake, expressed as a cyclic stress ratio; (2) the 
stiffness of the soil, expressed as an overburden stress-corrected shear wave velocity; and (3) the 
resistance of the soil to liquefaction, expressed as a cyclic resistance ratio. Each parameter is 
discussed below. 


2.1 CYCLIC STRESS RATIO (CSR) 


The cyclic stress ratio, tay /o',, at a particular depth in a level soil deposit can be 
expressed as (Seed and Idriss, 1971): 


Q 


CSRE V0 165 [2a [2 )r (2.1) 
Gy ig Ces 


where 


is = the average equivalent uniform shear stress caused by the earthquake and is 
assumed to be 0.65 of the maximum induced stress, 
= the peak horizontal ground surface acceleration, 


Omax 

g = the acceleration of gravity, 

Ov = the total vertical (overburden) stress at the depth in question, 

o', =. the initial effective overburden stress at the same depth, and 

ra = a shear stress reduction coefficient to adjust for the flexibility of the soil 
profile. 


2.1.1 Peak Horizontal Ground Surface Acceleration 


Peak horizontal ground surface acceleration is a characteristic of the ground shaking 
intensity, and is defined as the peak value in a horizontal ground acceleration record that would 
occur at the site without the influence of excess pore-water pressures or liquefaction that might 
develop (Youd et al., 2001). 


Peak accelerations are commonly estimated using empirical attenuation relationships of 
QAmax aS a function of earthquake magnitude, distance from the energy source or surface project 
of the fault rupture, and local site conditions. Since many published attenuation relationships are 
based on both peak values obtained from ground motion records for the two horizontal 
directions (sometimes referred to as the randomly orented horizontal component), the 
geometric mean (square root of the product) of the two peak values is used. According to 
Youd et al. (2001), use of the geometric mean is consistent with the derivation of the SPT-based 
procedure and is preferred for use in engineering practice. However, use of the larger of the 
two horizontal peak accelerations would be conservative and is allowable. 


- Regional or national seismic hazard maps (http://geohazards.cr.usgs.gov/eq/; Frankel et 
al., 2000) are also often used to estimate peak accelerations. If peak acceleration is estimated 
from a map, the magnitude and distance information should be obtained from the deaggregated 
matrices used to develop the map. The value of Qmax selected will depend on the target level of 
risk and compatibility of site conditions. For site conditions not compatible with available 
probabilistic maps or attenuation relationships, the value of a,,.. may be corrected based on 
dynamic site response analyses or site class coefficients given in the latest building codes. 


2.1.2 Total and Effective Overburden Stresses 


Required in the calculation of g,, and g', are densities of the various soil layers, as well 
as characteristics of the ground water. For non-critical projects involving hard-to-sample soils 
below the ground water table, densities are often estimated from typical values for soils with 
similar grain size and penetration or velocity characteristics. Fortunately, CSR is not very 
sensitive to density, and reasonable estimates of density yield reasonable results. 


The values of o', and CSR are sensitive to the ground water table depth. Other ground 
water characteristics that may be significant to liquefaction evaluations include seasonal and 
long-term water level variations, depth of and pressure in artesian zones, and whether the water 
table is perched or normal. 


2.1.3 Stress Reduction Coefficient 


Equation 2.1 is based on Newton’s second where force is equal to mass times 
acceleration. The coefficient rz is added because the soil column behaves as a deformable body 
rather than a rigid body. 


2.1.3.1 Relationship by Seed and Idriss (1971)—Values of rg are commonly 
estimated from the chart by Seed and Idriss (1971) shown in Fig. 2.1. This chart was 
determined analytically using a variety of earthquake motions and soil conditions. Average rq 
values given in the chart can be estimated using the following functions (Liao and Whitman, 
1986; Robertson and Wride, 1997): 


rg = 1.0 - 0.00765 z forz<9.15m (2.2a) 

rag= 1.174 - 0.0267 z for9.15m<z<23m (2.2b) 

ra = 0.744 - 0.008 z for 23 m<z<30m (2.2c) 
where 

Zz = the depth below the ground surface in meters. 


Figure 2.1 shows the average rz values approximated by Eq. (2.2). These average rz values 
were suggested by the 1996 NCEER Workshop (Youd et al., 1997; 2001) for non-critical 
projects. 


2.1.3.2 Revised Relationship Proposed by Idriss (1998; 1999)—Figure 2.2 presents 
revised average values of rz proposed by Idriss (1998; 1999) for various earthquake magnitudes. 
The plotted curves are averages of many individual curves derived analytically by Golesorkhi 
(1989) under the supervision of the late Prof. H. B. Seed. They are defined by the following 
relationship (after Idriss, 1998; modified for depth in meters): 


In(ra) = a(z) + (z) M, (233) 
where 

a(z) = -1.012 - 1.126 sin( +5.133), and (2.4) 

fz) = 0.106 + 0.118 sin +5.142) (2.5) 
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Fig. 2.1 - Relationship Between Stress Reduction Coefficient and Depth Developed by Seed 
and Idriss (1971) with Approximate Average Value Lines from Eq. 2.2. (after 


Youd et al., 1997) 
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Fig. 2.2 - Relationship Between Average Stress Reduction Coefficient and Depth Proposed 
by Idriss (1998; 1999) with Average of Range Determined by Seed and Idriss 
(1971). 


As shown in Fig. 2.2, the curve defined by Eq. 2.3 for M, = 7.5 is almost identical to the 
average of the range published by Seed and Idriss (1971). 


The scatter in the individual curves used to determine the average curves shown in Fig. 
2.2, as well as Fig. 2.1, is rather large. For example, coefficients determined for a 30 m thick, 
loose sand deposit and magnitude 5.5 earthquakes exhibit standard deviations of about 0.1 at a 
depth of 5 m and 0.15 at a depth of 10 m. These standard deviation values would be larger if 
soil deposits of various thicknesses and stiffnesses are considered. Thus, as an alternative 
approach, the variation of rz with depth may be calculated analytically using site-specific layer 
thicknesses and stiffnesses. 


2.2 STRESS-CORRECTED SHEAR WAVE VELOCITY 


The in situ Vs can be measured by several seismic tests including crosshole, downhole, 
seismic cone penetrometer (SCPT), suspension logger, and Spectral-Analysis-of-Surface-Waves 
(SASW). Recent reviews of these test methods are given in Woods (1994), Kramer (1996), and 
Ishihara (1996). ASTM D-4428M-91 provides a standard test method for crosshole seismic 
testing. Standard test methods do not exist for the other seismic tests. Primary advantages and 
disadvantages of the various in situ Vs test methods are presented in Table 2.1. The accuracy of 
each test method can be sensitive to equipment and procedural details, soil conditions, and 
interpretation techniques. 


One important factor influencing Vs is the state of stress in soil (Hardin and Dmevich, 
1972; Seed et al., 1986). Laboratory test results (Roesler, 1979; Stokoe et al., 1985; Belloti et 
al., 1996) show that the velocity of a propagating shear wave depends equally on principal 
stresses in the direction of wave propagation and the direction of particle motion. Thus, Vs 
measurements made with wave propagation or particle motion in the vertical direction can be 
related by the following empirical relationship: 


Vs=A (o')" (o',)” (2.6) 
where 

A = a parameter that depends on the soil structure, 

o', = the initial effective horizontal stress at the depth in question, and 

m = a Stress exponent with a value of about 0.125. 


10 


‘spoyjoul []e@ 10 Sutseds 4so} UO spuodap uonnjosoy, 


‘spoyjour [je J0J sonbruyos) uoneyosdioqut poos pur “sjreyop jeinpoooid pure yuouldinbo poos uo spuadap Aq1yenb poo, 


S]IOs 


ayenouod 0} INOWIp QU09 YIM a]1yOId 

dunso} pue seose parade] poylejyop ‘ou0o 
spose o31e] Surusosos | odsey Suideumt 1ydesd OMUSIOS WO Poure}qo poyjow 
IO} PojINs-[]9\4. -OWO} JOJ POJINS-[[9\\ osje eyep UONeqUed 4S9} 9]QRI]II SOW] 


(ZH 0002 
-ZH 00s) Aousnboyy 


389} YsTY ‘sTIOs 
YOs pue sojoy posed ul 
yiom jou Aew ‘pormbes 


soy por[y-PmMy 


poJoA0001 sojdures ou 
‘s1oAe] INS 0} JudONfpe 
jewojeul yos pue 

sXe] UN JO UOTNNIOSEI 
1ood ‘pounsse 
BuLsAP] [eJUOZLOP] 


poyjou 
1S9} SAISUSdXO SOU 
‘yydop 3s0} 78 [BLAVeUr 
Jays sasuas ‘swojqoid 
UONIeIJOI 9{QISSOd 


undap YyIM soseosOUl 
yjed foaey oar ‘SIdAe] 
MOTTEYS ILM suso{qosd 
UONORIJAI JIQISSOd 


posoaoool sojdures ou 
‘undop wpim surseoioul 
AJLOOJ9A JOF SYIOM 
ATuo ‘jso} UONeIJOI UY 


= 


UONOIIIP [ROIOA 

ul uonededoid sABM JO 
uonow ojo1ed 

Jo juauodusoo soley 


QoINOs ureDq JBSYs 
puke a]0y210q U90A139q 
QOURISIP URY} Ja}BOId 
yidop 3s9} ym ‘Sox 


UONOBIJOI 
JOJ SONOS [POA pure 
UOT}OaYJaI OJ 20N0s 
jeOZLOY WIM ‘sax, 


SOARM IBOYs pozuLjod 
AT[BOTOA YEA “SOX 


90INnOS 
JROHOA YIM ‘SOX 


odBI9AB [eqO|s poos 
sopraoid ‘uyydap yy 
SOSBOIDOP IIB} 0} POOH 


Aipiqeuea ssougns 
[IOs JO UOTNTOSaYy 


QOBpINS 9} 0} SSOJD yydop WIM yidop 
1ood ‘yjdap 38 pooy | sasealoop *11eJ 0} POOH YIM JUR}SUOD “Pood 


porinbor (sdulpunos 
I QJOUI 10 Z JO) SoTOYSI0G JO ‘ON 


I9}9WI0J}QUIg BUOD 9]0YSsOlD) | 
OMUSI9g 7 2]0YUMOK] 


‘jUoWUSsassy UONORJaNbry JO} SpoyIaP ISP 54 MIS UT SNOLeA Jo sosejueapesiq pue sosejueapy Jo uosuedwoy - [7 IqeL 


QSBIOAR [BQO[S BSIBOD 
sopraoid ‘100d 0} 11e.J 


yeauue 
dABM IBOYS YSINsUNsIp 
0} J[NOIIp usyo ‘11e.] 


Aniqeyeadas 


enbiuyos} uonejoidiojur 
pue joNuod AVITENd 


xa]dw09 *11ey 0} pooh 


UOIPOROY 
fuonoayoy soejing 


/ SOARA\-BdRJING 
jo-sisAyeuy-[eroeds 


uorsuodsng 


poyjepy JoWoInsesjy 


11 


Following the traditional procedures for correcting SPT blow count and CPT tip 


resistance, one can correct Vs to a reference overburden stress by (Sykora, 1987b; Robertson et 
al., 1992): 


0.25 
Vs, = Vs Cys =Vs (2a (277) 
o'y 
where 
Vs; = the overburden stress-corrected shear wave velocity, 
Cys = a factor to correct measured shear wave velocity for overburden pressure; 
Pa = areference stress, 100 kPa or approximately atmospheric pressure, and 
o', =. the initial effective overburden stress in kPa. 


A maximum Cys value of 1.4 is generally applied to Vs data at shallow depths, similar to the SPT 
and CPT procedures. In using Eq. (2.7), it is implicitly assumed that the initial effective 
horizontal stress, o',, 1s a constant factor of the initial effective overburden stress (because both 
o', and go', affect Vs as shown in Eq. (2.6)). This factor, generally referred to as K%, is 
assumed to be approximately 0.5 at natural, level-ground sites where liquefaction has occurred 
or is likely to occur. Also, in applying Eq. (2.7), it is implicitly assumed that Vs is measured with 


both the directions of particle motion and wave propagation polarized along principal stress 
directions and one of those directions is vertical. 


Since the direction of wave propagation and the direction of particle motion is different 
with respect to the stress in the soil for each in situ seismic test method, some variations 
between measured Vs is expected. These variations are minimized by performing the tests with 
at least a major component of wave propagation or particle motion in the vertical direction. To 
have a major component of wave propagation or particle motion in the vertical direction, 
crosshole tests are conducted with particle motion in the vertical direction, downhole and 
seismic cone tests are conducted at depths greater than the distance between the shear beam 
source and the borehole or cone sounding such that wave propagation is in the vertical direction, 
and SASW tests are conducted with a vertical source. 


In soils above the ground water table, particularly silty soils, negative pore pressures 
increase the effective state of stress and, hence, the value of Vs; measured in seismic tests. This 


effect should be considered in the estimation of o', for correcting Vs to Vs;, and for computing 
~ CSR using Eq. (2.1). 
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2.3 CYCLIC RESISTANCE RATIO (CRR) 


The value of CSR separating liquefaction and non-liquefaction occurrences for a given 
Vs}, or corrected penetration resistance, is called the cyclic resistance ratio, CRR. Figure 2.3 
presents the CRR-Vs; curves developed by Andrus et al. (1999) for magnitude 7.5 earthquakes 
and uncemented, Holocene-age soils. The curves are dashed above CRR of about 0.35 to 
indicate that they are based on limited field performance data, as discussed in Appendix F. The 
curves do not extend much below 100 m/s, since there are no field data to support extending 
them to the origin. They are defined by: 


2 
CRR = MSF 0.o22{ Kshs) “3 2 -=| (2.8) 
100 Vsig Kel st is! 
where 
MSF = the magnitude scaling factor to account for the effect of earthquake 
magnitude, 

Vs, = the limiting upper value of Vs, for cyclic liquefaction occurrence, and 
K, = a factor to correct for high Vs; values caused by cementation and aging. 


The first (or squared) term in Eq. (2.8) is based on a relationship between CRR and Vs, for 
constant average cyclic strain derived by R. Dobry. The second term is a hyperbola with small 
value at low values of Vs;, and a very large value as Vs; approaches )’5,. 


2.3.1 Magnitude Scaling Factor 


The magnitude scaling factor is traditionally applied to CRA, rather than the cyclic 
loading parameter CSR, and equals 1 for earthquakes with a magnitude of 7.5. For magnitudes 
other than 7.5, Table 2.2 presents scaling factors developed by various investigators. These 
magnitude scaling factors were derived from laboratory test results and representative cycles of 
loading (Seed and Idriss, 1982; Idriss, personal communication to T. L. Youd, 1995; Idriss, 
1998; Idriss, 1999), correlations of field performance data and blow count measurements 
(Ambrasey, 1988; Youd and Noble, 1997), estimates of seismic energy for laboratory and field 
data (Arango, 1996), and correlations of field performance data and in situ Vs; measurements 
(Andrus and Stokoe, 1997). Figure 2.4 shows a plot of the various magnitude scaling factors 
along with the range recommended by the NCEER Workshop (Youd et al., 1997; 2001). 
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Fig. 2.3 - Curves Recommended for Calculation of CRR from Vs, (Andrus et al., 1999). 


Table 2.2 - Magnitude Scaling Factors Obtained by Various Investigators. (modified from 
Youd and Noble, 1997) 


Magnitude Scaling Factor (MSF) 


Ambraseys Youd & Nobie Arango | Andrus 
(1988) (1997) (1996)** & 

P,, % Stokoe 

<20 <30 <50 (1997) 


*Extrapolated from scaling factors for M, = 6, 6.5, 7, and 7.5 using MSF = (M,/7.5)>3. 
**Based on equivalent uniform number of cycles and consideration of distant liquefaction sites (Column 10), 


and energy principles (Column 11). 


Although the 1996 NCEER Workshop (Youd et al., 1997) recommended a range of 
magnitude scaling factors for engineering practice, a consensus has not yet been reached by the 
workshop participants. At the August 1998 MCEER Workshop, a revised set of magnitude 
scaling factors and stress reduction coefficients (see Section 2.1.3.2) were proposed by I. M. 
Idriss. Liao and Lum (1998) present results of a statistical analysis supporting the original Seed 
and Idriss (1982) factors. The magnitude scaling factors recommended by the 1996 NCEER 
Workshop and the revised factors proposed by Idriss (1999) are discussed below. 


2.3.1.1 Factors Recommended by 1996 NCEER Workshop—tThe magnitude scaling 
factors recommended by the 1996 NCEER Workshop (Youd et al., 1997) can be represented 


by: 


Mw) 
5. el bapa 2.9 
sane ( 1) ) ae 
where 
M, = moment magnitude, and 


an exponent. 


= 
Il 
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Magnitude Scaling Factor, MSF 


Range of Recommended 
MSF from the 1996 
NCEER Workshop 


Seed & Idriss (1982) 
Idriss (1995) 
Idriss (1998) 
Idriss (1999) 
Ambraseys (1988) 
Arango (1996) 
Arango (1996) 
Andrus & Stokoe (1997) 
Youd & Noble (1997) 

Pr < 20 % 

Pa < 32 

P. < 50 


Range of Recommended 
MSF from the 1996 
NCEER Workshop 


Earthquake Magnitude, My, 


Fig. 2.4 - Magnitude Scaling Factors Derived by Various Investigators with Range 


Recommended by the 1996 NCEER Workshop (after Youd et al., 1997; 
Youd and Noble, 1997). 


Moment magnitude is the scale most commonly used for engineering applications, and is 
preferred for liquefaction resistance calculations (Youd et al., 1997). When only other 


magnitude scales are available, they can be converted to M,, using the relationship of Heaton et 
al. (1982) shown in Fig. 2.5. 


The lower bound for the range of magnitude scaling factors recommended by the 1996 
NCEER Workshop is defined with n = -2.56 (Idriss, personal communication to T. L. Youd, 
1995) for earthquakes with magnitudes < 7.5. The upper bound of the recommended range is 
defined with n = -3.3 (Andrus and Stokoe, 1997) for earthquakes with magnitudes < 7.5. For 
earthquakes with magnitudes > 7.5, the recommended factors are defined with n = -2.56. 
Magnitude scaling factors defined by Eq. (2.9) and average rz values originally proposed by 
Seed and Idriss (1971) should be used together when applying Eqs. (2.1) and (2.8). 


2.3.1.2 Revised Factors Proposed by Idriss (1999)—The magnitude scaling factors 
proposed by Idriss (1999) are derived using laboratory data from Yoshimi et al. (1984) and a 
revised relationship between representative cycles of loading and earthquake magnitude. They 
are defined by the following equation: 


MSF = 6.9e9{ = =| — 0.06 for M, > 5.2 (2.10a) 


MSF = 1.82 for M, < 5.2 (2.10b) 


where exp is the constant e raised to the power of the number given in the parentheses. 
Magnitude scaling factors defined by Eq. (2.10) and revised 7g proposed by Idriss (1999) should 
be used together when applying Eqs. (2.1) and (2.8). 


2.3.1.3 Recommended Magnitude Scaling Factors—There is little difference in using 
magnitude scaling factors and rz values recommended by the 1996 NCEER Workshop (Youd et 
al., 1997) and those proposed by Idriss (1999) for magnitudes near 7.5 and depths less than 11 
m (see Appendix F, Sections F.2.2 and F.2.3). At magnitudes near 5.5, the difference is 
significant. The lower bound for the range of magnitude scaling factors defined by Eq. (2.9) 
with m = -2.56 is recommended in these guidelines because it provides more conservative 
assessment than with n = -3.3 for magnitudes less than 7.5. While the magnitude scaling factors 
defined by Eq. (2.9) with 7 = -2.56 are less conservative than the factors proposed by Idriss 
(1999) for magnitudes less than 7.5, the findings of Ambrasey (1988), I. M. Idnss (personal 
communication to T. L. Youd, 1995), Arango (1996), Youd and Noble (1997), and Andmus and 
Stokoe (1997; as indicated by the very conservative CRR-Vs; curves shown in Figs. F.13 and 
F.14) support their use. 
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Fig. 2.5 - Relationship Between Moment Magnitude and Various Magnitude Scales 
(after Heaton et al., 1982). 
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2.3.2 Limiting Upper Value of Vs; in Sandy Soils 


The assumption of a limiting upper value of Vs; is equivalent to the assumption 
commonly made in the SPT- and CPT-based procedures dealing with clean sands, where 
liquefaction is considered not possible above a corrected blow count of about 30 (Seed et al., 
1985) and a corrected tip resistance of about 160 (Robertson and Wride, 1998). Upper limits 
for Vs; and penetration resistance are explained by the tendency of dense soils to exhibit dilative 
behavior at large strains, causing negative pore-water pressures. While it is possible in a dense 
soil to generate pore-water pressures close to the confining stress if large cyclic strains or many 
cycles are applied, the amount of water expelled during reconsolidation is dramatically less for 
dense soils than for loose soils. As explained by Dobry (1989), in dense soils, settlement is 
insignificant and no sand boils or failure take place because of the small amount of water 
expelled. This is important because the definition of liquefaction used to classify the field 
behavior here, as well as in the penetration-based procedures, is based on surface manifestations 
such as boils and ground cracks. 


The case history data above a CSR value of about 0.35 are limited, as discussed in 
Appendix F. Thus, current estimates of ’s, rely, in part, on penetration-Vs correlations and, in 
part on the case histories. Values of 5, can be estimated from: 


Vs, = 215 m/s for sands with FC <5 % (2.11a) 

Vis) = 215 -0.5(FC-5) m/s for sands with 5% < FC <35% (2.11b) 

V5; = 200 m/s for sands and silts with FC > 35 % (2.11c) 
where 

FC = average fines content in percent by mass. 


Equations (2.8) and (2.11a) provide a CRR value of about 0.6 at Vs; = 210 m/s. A Vs; value of 
210 m/s is considered equivalent to a corrected blow count of 30 in sands with FC = 5 %, based 
on penetration-Vs correlations. 


2.3.3 Limiting Upper Value of Vs; in Gravelly Soils 


Although the )’5; values given in Eq. (2.11) were determined for sandy soils, the case 
history data indicate that these limits also represent reasonable limits for gravelly soils divided 
into the same categories based on fines content (see Fig. F.7). This might be considered rather 
surprising based on the penetration-Vs correlations presented in the literature for gravelly soils. 
For instance, the correlation by Ohta and Goto (1978) suggested a Vs; value of 227 m/s for 
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Holocene gravels at an equivalent (N1)60 of 30. Similarly, the correlation by Rollins et al. (1998) 
provided a best-fit Vs; value of 232 m/s for Holocene gravels. On the other hand, all the 
liquefaction case history data exhibit Vs; values of about 200 m/s or less, suggesting that 230 
m/s may be inappropriately high. To investigate further the value of Vs; in gravelly soils, 
laboratory studies involving Vs; measurements in gravelly soils were reviewed. Kokusho et al. 
(1995) clearly showed that the shear wave velocity (or stiffness) of gravelly soils varies greatly 
and is highly dependent on the particle gradation. Weston (1996) showed similar results for 
coarse sands with gravels. In both cases, the results show that increasing the uniformity 
coefficient can significantly increase the shear wave velocity in medium dense to dense gravels. 
On the other hand, very loose gravelly soils, even well-graded gravels, can exhibit shear wave 
velocities similar to those of loose sands (Kokusho et al., 1995). The case history data 
presented in Fig. F.7 support the premise that gravelly soils that are loose enough to exhibit 
significant liquefaction effects (boils, ground cracks, etc.) have shear wave velocities similar to 
loose sands. Hence, the boundaries developed for sandy soils are recommended as preliminary 
boundaries for gravelly soils. However, additional work is clearly needed to understand the 
relationship between Vs, and liquefaction resistance of gravels. 


2.3.4 Cementation and Aging Correction Factor 


The recommended CRR-Vs; curves shown in Fig. 2.3 are limited to the characteristics of 
the database used to develop them. The database consists of relatively level ground sites with 
the following general characteristics: (1) uncemented soils of Holocene age; (2) average depths 
less than about 10 m; (3) ground water table depths between 0.5 m and 6 m, and (4) all Vs 
measurements are from below the water table. Correction factors may be used to extend the 
curves to site conditions different from the database. 


The correction factor K, is 1 for areas of uncemented, Holocene-age soils. For 
Pleistocene-age soils (>10 000 years), average estimates of K, range from 0.6 to 0.8 based on 
the penetration-Vs; correlations by Rollins et al. (1998a) and Ohta and Goto (1978), 
respectively. Figures 2.6 and 2.7 illustrate two methods for estimating the value of K, using 
SPT and CPT test results, respectively. Shown in the figures are correlations for clean sands 
and silty soils implied by the CRR-Vs, relationship defined by Eq. (2.8) and 1996 NCEER 
Workshop recommended CRR-penetration relationships (Youd et al., 1997). In the example, 
the measured values of Vs, (Ni)6o, Gcin, and fines content are 220 m/s, 8, 55, and 10 %, 
respectively. The Vs;-penetration correlations in Figs. 2.6 and 2.7 suggest a K, value between 
0.71 and 0.75, respectively, for these conditions. The K, value is assumed to be the ratio of the 


predicted value of Vs;, based on the corrected penetration resistance and fines content, to the 
measured value of V5;. 
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Fig. 2.6 - Suggested Method for Determining the Correction Factor K; from (Nj)69, Vs;, and 


Fines Content at a Weakly Cemented Soil Site. 
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Fig. 2.7 - Suggested Method for Determining the Correction Factor K, from q,j;y, Vs), and Fines 
Content at a Weakly Cemented Soil Site. 
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The method for estimating K, described above assumes that the strain level induced 
during penetration testing is the same strain level causing liquefaction, which may not be true 
because pore-water pressure buildup to liquefaction can occur at medium strains in several 
loading cycles (Dobry et al., 1982; Seed et al, 1983). The method also assumes that 
liquefaction potential, blow count, and cone penetration resistance are not affected by 
cementation, which may not be a reasonable assumption. Hence, this suggested method should 
be used cautiously and with engineering judgment. 


ue 


2.4 FACTOR OF SAFETY 


A common way to quantify the potential for liquefaction is in terms of a factor of safety. 
The factor of safety, F's, against liquefaction can be defined by: 


_ CRR 
CSR 


F's (2212) 


By convention, liquefaction is predicted to occur when Fs < 1. When Fs > 1, liquefaction is 
predicted not to occur. 


As is the case with the SPT- and CPT-based charts, it is possible that liquefaction could 
occur outside the region of predicted liquefaction shown in Fig. 2.3. Consequently, the Building 
Seismic Safety Council (1997, page 158) suggests a factor of safety of 1.2 to 1.5 is appropriate 
when applying the Seed-Idriss simplified procedure in engineering design. The acceptable value 
of F's for a particular site will depend on several factors, including the type and importance of 
structure and the potential for ground deformation. Based on Vs-SPT blow count correlations 
(see Section F.4) and probability studies (see Section G.2), the recommended Vs-based 
procedure is as conservative as the SPT-based procedure outlined by Seed et al. (1985) and 
updated by the NCEER Workshop (Youd et al., 2001). Thus, the same range of factor of safety 
is recommended for the Vs-based method. 


2.5 PROBABILITY-BASED EVALUATION 


Probability of liquefaction, P,, is required information for making risk-based design 
decisions. As discussed in Appendix G, the relationship between P; and Fs for the deterministic 
procedure described above can be expressed as (Juang et al., 2002; modified from Juang et al., 
2001a): 


Pe 


j CAEN 


In Eq. (2.13), a Fs value of 1 corresponds to points on the deterministic curves shown in Fig. 
2.3. Thus, on average, the deterministic curves are characterized by a P, value of 26 %. This 
average P, value is similar to probability estimates determined for the SPT-based curves (Liao et 
al., 1988; Youd and Noble, 1997; Juang et al., 2000a). The relationship defined by Eq. (2.13) is 
plotted in Fig. 2.8, and provides the link between the probabilistic and deterministic methods. 
By combining Eqs. (2.8), (2.12) and (2.13), one can obtain a family of P; curves for risk-based 
design. The family of P, curves for magnitude 7.5 earthquakes and soils with FC < 5 % is 
presented in Fig. 2.9. 


It is important to note that Figs. 2.8 and 2.9 are developed assuming F's to be a fixed 
variable, and possible variations in CRR and CSR are not considered directly. Previous studies 
by Juang et al. (2000b) and Chen and Juang (2000) concluded that practically the same P,-F's 
relationship would be obtained even if the uncertainties in CSR and CRR were incorporated in 
the formulation for P;. Thus, in general, if the variations of CRR and CSR are not too great, the 
figures can be used directly without considering the variations (Juang et al., 2001b). 


2.6 SUMMARY 


. In this chapter, guidelines are presented for evaluating liquefaction resistance through Vs 
measurements using the procedure outlined in Andrus and Stokoe (2000). The procedure can 
be summarized in the following ten steps: 


1. From available subsurface data, develop detailed profiles of Vs, soil type, fines 
content and, if possible, soil density and penetration resistance. Identify the 
depth of the ground water table, noting any seasonal fluctuations and artesian 
pressures. 


2. Calculate the values of g, and o', for each measurement depth at which 
seismic testing has been performed. 


3. Correct the Vs measurements to the reference overburden stress of 100 kPa 


using Eq. (2.7). The correction factor Cys is limited to a maximum value of 
1.4 at shallow depths. 
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Relationship based on: 
Andrus et al. (1999) procedure 


and case histories 
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probability density functions of 
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Fig. 2.8 - Suggested Relationship for Selecting F's Based on Probability of Liquefaction and the 
Recommended CRR-Vs; Curves. (Juang et al., 2002) 
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4. Determine the value of )’s, for each measurement depth using Eq. (2.11) 
which is recommended for sandy as well as gravelly soils. If the fines content 
is unknown, assume 215 m/s for V5). 


5. Determine the value of K,. K, can be assumed equal to 1, if the soil to be 
evaluated is uncemented and less than 10 000 years old. If the soil conditions 
are unknown and penetration data are not available, assume 0.6 for Kz. 


6. Determine the design earthquake magnitude and expected value of Qmax. 


7. Calculate CSR for each measurement depth below the water table using Eq. 
(2.1). The value of rz can be estimated from the average curve originally 
proposed by Seed and Idriss (1971). 


8. Calculate CRR for each value of Vs; using Eqs. (2.8) and (2.9) with n = -2.56. 
It is important to note that Eq. (2.8) is for extreme behavior where boils and 
ground cracks occur. 


9. Calculate the value of Fs for each value of Vs, using Eq. (2.12). By 
convention, liquefaction is predicted to occur when Fs < 1, and not to occur 
when F's > 1. 


10. Plot the values of Vs;, CSR, CRR and Fs to visually note how they vary with 
depth, and how many points fall in the regions of liquefaction and no 
liquefaction. 


The deterministic Vs-based procedure outlined above is characterized with an average 
probability of liquefaction of 26 %. In other words, a soil with a calculated F's = 1 has a 26 % 
chance of liquefaction occurrence based on the case histories analyzed in this study. As 
mentioned previously, the Vs-based procedure is as conservative as the SPT-based procedure by 
Seed et al. (1985). A factor of safety of 1.2 to 1.5, as suggested by the Building Seismic Safety 
Council (1997, page 158), is considered appropriate for design of typical buildings using the 
SPT- and Vs-based procedures. This range of factor of safety corresponds to a probability of 
liquefaction of about 8 % (for Fs = 1.5) to 16 % (for Fs = 1.2). The acceptable value of F's for a 
particular site will depend on several factors, including the type and importance of structure and 
the potential for ground deformation. For critical structures, a smaller probability of liquefaction 
might be required. Equation (2.13) provides an important link between Fs and P,, and is 
suggested for probability liquefaction evaluations. 
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CHAPTER 3 


APPLICATION OF THE LIQUEFACTION EVALUATION PROCEDURE 


To illustrate the application of the liquefaction evaluation procedure described in Chapter 
2, two sites shaken by the 1989 Loma Prieta, California, earthquake (M4, = 7) are considered 


below. The two sites are Treasure Island Fire Station and Marina District Winfield Scott 
School. 


3.1 TREASURE ISLAND FIRE STATION 


Treasure Island is a man-made island located in the San Francisco Bay along the Bay 
Bridge between the cities of San Francisco and Oakland. It was constructed in 1936-37 by 
hydraulic filling behind a perimeter rock dike. The perimeter dike served to contain the 
hydraulic fill and was raised in sections over the previously placed fill. In 1991, Treasure Island 
was selected as a national geotechnical experimentation site. Much of the work to date at the 
Treasure Island national geotechnical experimentation site centers arounds a ground response 


experiment (de Alba and Faris, 1996) with six accelerometers and eight piezometers operating at 
various elevations near the fire station. 


Extensive field tests have been conducted near the Treasure Island fire station to 
characterize ground conditions. Figures 3.1 and 3.2 present Vs and general soil profiles for the 
site. The Vs profile shown in Fig. 3.1(a) is from Fuhriman (1993), and was determined by 
crosshole testing. The Vs profile shown in Fig. 3.2(a) is based on unpublished SASW test results 
by The University of Texas at Austin in 1992. From the description by de Alba et al. (1994), the 
upper 4.5 m of soil consists of silty sand fill, possibly formed by dumping. Between depths of 
4.5 m and 12.2 m, the soil consists of silty sand to clayey sand, formed by hydraulic filling. 
Beneath the hydraulic fill are natural clayey soils. The ground water table lies near the ground 
surface at a depth of 1.4 m. 


During the 1989 Loma Prieta earthquake, a seismograph station at the fire station 
recorded ground surface accelerations. Unlike recordings at other seismograph stations located 
on soft-soils in the Bay area, there is a sudden drop in the recorded acceleration at about 15 
seconds and small motion afterward (Idriss, 1990). De Alba et al. (1994) attribute this behavior 
to liquefaction of an underlying sand layer, although no sand boils or ground cracks occurred at 
the site. The nearest liquefaction effect observed is a sand boil located 100 m from the site 
(Geometric Consultants, 1990; Bennett, 1994; Power et al., 1998). 
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Figure 3.1 presents the liquefaction evaluation for the crosshole test array B1-B4 and the 
1989 Loma Prieta earthquake. Values of Vs; and CSR shown in Figs. 3.1(a) and 3.1(d), 
respectively, are calculated assuming densities of 1.76 Mg/m’ above the water table and 1.92 
Mg/m®* below the water table. Based on peak values of 0.16 g and 0.11 g recorded in two 
horizontal directions at the fire station dunng the 1989 earthquake (Brady and Shakal, 1994), a 
geometric mean value of 0.13 g is used to calculate CSR. Stress reduction coefficients are 
estimated using the average curve by Seed and Idriss (1971) shown in Fig. 2.1. 


For the crosshole measurement at a depth of 4.6 m, values of CSR and Vs; are calculated 
as follows: 


CSR =0.65| 22% ]{ 2 |,, = 0.65] 2138 |( 849) 97 = 0.131 (3.1) 
g Jlo g )\527 


and 


0.25 
pe 100 
Vs; = (2 7 = 134( 200)" = 158 m/s (3.2) 


Assuming an average fines content of 24 %, from Fig. 3.1(c), and a K, value of 1, the values of 
Vs}, CRR, and Fs are calculated by: 


Vx = 215 - 0.5(FC-5) = 215 - 0.5(24-5) = 206 m/s (3.3) 
and 
aN 1 1 
CRR = of X21) {4-5 MSF (3.4) 
100 Va-KVs Vs 
2.56 
- {o.022( 155) 12a (=) 
100 206-158 206)( \7.5 
= 0.119 
and 
CRR 0.119 
= —— =—— = 09] 
ty GSR 31 oy 


Since the value of F's is less than 1, liquefaction is predicted at this depth. 
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Values of Fs shown in Fig. 3.1(e) are less than 1 for the depths of 4 m to about 9 m. 
Between the depths of 4 m and 7 m, the sand contains non-plastic fines and is considered 
liquefiable. Between the depths of 7 m and 9 m, the soil exhibits plastic characteristics and may 
be non-liquefiable by the simple clay criteria (see Section 1.1). Thus, the layer most likely to 
liquefy, or the critical layer, lies between the depths of 4 m and 7 m. 


Figure 3.2 presents the liquefaction evaluation for the SASW test array. Locations of Vs 
measurements for the SASW test array are assumed at the center of the layer used in forward 
modeling of surface wave measurements. Values of Fs shown in Fig. 3.2(e) are less than 1 
between the depths of about 3.5 m and 11 m. The lowest values of Fs in the non-plastic soil is 
0.75 at a depth of 5.3 m. This F's value is similar to the lowest F's value of 0.77 determined from 
crosshole measurements in the critical layer. 


Figures 3.3 and 3.4 present the liquefaction evaluations directly on the recommended 
liquefaction assessment chart for the crosshole test array and SASW test array, respectively. 
Plots of this type are particularly useful in comparing the range and distribution of Vs, and CSR 
values with the case histories used to develop the assessment chart. Based on Fig. F.15, the 
assessment chart is well supported within the range of the Treasure Island Fire Station data. 


Although no sand boils or ground cracks occurred at the fire station during the 1989 
earthquake, the prediction of liquefaction agrees with the conclusion stated above that 
liquefaction of an underlying sand cause the sudden drop in the acceleration time histories 
recorded at this site (de Alba et al., 1994). A similar sudden drop in the strong ground motion 
recordings occurred at the Port Island Downhole Array site in Kobe, Japan, during the 1995 
Hyogo-ken Nanbu earthquake (Aguirre and Irikura, 1997), where liquefaction and sand boils did 
occur. It is possible that the 4 m thick layer capping the site, predicted not to liquefy, as shown 


in Figs. 3.1(e) and 3.2(e), prevented the formation of sand boils at the ground surface (Ishihara, 
1985). 


3.2 MARINA DISTRICT WINFIELD SCOTT SCHOOL 


Kayen et al. (1990) conducted downhole seismic tests at the Winfield Scott School in the 
Marina District of San Francisco. Figures 3.5(a) and 3.5(b) present the Vs and soil profiles for 
the site. The Vs profile was originally determined based on best-fit line segments through travel 
time measurements plotted versus depth. However, the layering assumed in the best-fit segment 
method did not seem appropriate for the fill. Figures 3.5(a) presents the Vs profile for the site 
determined using the pseudo-interval method (see Appendix E). Figure 3.5(c) presents a profile 
of fines content that are based on information provided by Kayen et al. (1990). The upper 7.6 m 
of soil at the site consists of sand with 1 % to 8 % fines. The ground water table lies at a depth 
of 2.7 m. 
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Chart Based on: AA ee7, 
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Content (%) 


Cyclic Stress or Resistance Ratio, CSR or CRR 
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Overburden Stress-Corrected Shear Wave 
Velocity, Vs7, m/s 


Fig. 3.3 - Liquefaction Assessment Chart for Magnitude 7 Earthquakes with Data for 
the 1989 Loma Prieta Earthquake and the Treasure Island Fire Station Site, 
Crosshole Test Array B1-B4 (Depths of 1.5 m to 14 m). 
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Chart based on: My = 7 

Average values of eae 
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Fig. 3.4 - Liquefaction Assessment Chart for Magnitude 7 Earthquakes with Data for 
the 1989 Loma Prieta Earthquake and the Treasure Island Fire Station Site, 
SASW Test Array (Depths of 2 m to 13 m). 
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Many structures, pavements, and public works near the school sustained heavy damage 
during the 1989 Loma Prieta earthquake (Kayen et al. 1990). This damage was due to 
liquefaction of the sand fill. From maps prepared by Pease and O’Rourke (1995), the site lies on 
the margin of the 1906 water front and artificial fill where about 40 mm of settlement occurred. 
Mapped sand boils and ground cracks lie just east of the site. Based on these observations, this 
site is classified as a liquefaction site during this earthquake. 


The Marina District and Treasure Island are located about 82 km from the 1989 surface 
fault rupture. Assuming a distance of 82 km from the fault rupture, the attenuation relationship 
by Idriss (1991) for 1989 strong ground motion records from soft-soil sites provides a median 
value of 0.16 g. This value is slightly higher than the geometric mean value of 0.13 g for the 
two peak horizontal accelerations recorded at Treasure Island fire station. Thus, a peak 
horizontal ground surface acceleration of 0.15 g, the average of these two estimates, is assumed 
in the analysis. 


Figure 3.5 presents the liquefaction evaluation for the Marina District School site and the 
1989 Loma Prieta earthquake. Values of Vs; and CSR are calculated assuming densities of 1.76 
Mg/m? above the water table and 1.92 Mg/m” below the water table. They are plotted in Fig. 
3.5(a) and 3.5(d) at the depths midway between receiver locations. Since the sand is uncemented 
and less than 10 000 years old, the value of K, is 1. Calculated values of F's are 0.42, 0.90, and 
0.51 at the depths of 3 m, 4 m, and 6.7 m within the sand fill. The silty clay layer beneath the 
sand fill is non-liquefiable by the simple clay criteria (see Section 1.1). Thus, having the lowest 
average value of Fs, the sand fill just below the water table between the depths of 2.7 m and 4.4 
m is identified as the critical layer that liquefied. A prediction of liquefaction agrees with the 
observed field behavior. 


Figure 3.6 present the liquefaction evaluation directly on the recommended liquefaction 
assessment chart for the Marina District. Based on Fig. F.15, the assessment chart is well 
supported within the range of the Marina District data. 
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Chart based on: Mirai 
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Overburden Stress-Corrected Shear Wave 
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Fig. 3.6 - Liquefaction Assessment Chart for Magnitude 7 Earthquakes with Data for 
the 1989 Loma Prieta Earthquake and the Marina District School Site (Depths 
of 3 mto7 m). 
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CHAPTER 4 


SUMMARY AND RECOMMENDATIONS 


4.1 SUMMARY 


Presented in this report are guidelines for evaluating the liquefaction resistance of soils 
through shear wave velocity, Vs, measurements. The guidelines are based on an earlier report 
entitled “Draft Guidelines for Evaluating Liquefaction Resistance Using Shear Wave Velocity 
Measurements and Simplified Procedures.” From comments received on the earlier report, the 
draft guidelines are updated in this report. The guidelines present a recommended procedure for 
evaluating soil liquefaction resistance and guidance for its use. 


The recommended procedure follows the general format of the simplified penetration- 
based procedure originally proposed by Seed and Idriss (1971). Cyclic stress ratios, CSR, are 
calculated using Eq. (2.1), with the average stress reduction coefficient estimated from Fig. 2.1. 
Shear wave velocity measurements are corrected for overburden stress using Eq. (2.7). Figure 
2.3 presents the recommended evaluation curves for uncemented, Holocene-age soils and 
magnitude 7.5 earthquakes. These curves are defined by Eq. (2.8) with MSF = 1, Vs, = 200 
m/s to 215 m/s (depending on fines content), and K, = 1. Equation (2.8) can be adjusted for 
other magnitude earthquakes using MSF values defined by Eq. (2.9) with n= -2.56. Corrections 


for cemented and aged soils are suggested in Section 2.3.4. A ten-step summary of the 
procedure is given in Section 2.6. 


The recommended liquefaction evaluation curves, defined by Eq. (2.8), are based on a 
modified relationship between overburden stress-corrected shear wave velocity, Vs;, and CSR 
for constant average cyclic shear strain suggested by R. Dobry. As discussed in Section 2.3 and 
Appendix F, the quadratic relationship proposed by Dobry is modified so that it is asymptotic to 
some limiting upper value of Vs;. This limit is related to the tendency of dense granular soils to 
exhibit dilative behavior at large strains, as well as the fact that dense soils expel dramatically 
less water during reconsolidation than loose soils. Liquefaction and non-liquefaction case 
histories from 26 earthquakes and more than 70 measurement sites in soils ranging from clean 
fine sand to sandy gravel with cobbles to profiles including silty clay layers are analyzed to 
determine the parameters of Eq. (2.8). Penetration-Vs correlations are also considered. The 
evaluation curves correctly bound over 95 % of the case histories where liquefaction occurred. 


an 


By constructing relationships between Vs, and penetration resistance from the recommended 
evaluation curves and plotting available in situ test data, it is shown that the Vs,-based 
evaluation curves are generally more conservative than the penetration-based evaluation curves. 
From logistic regression and Bayesian interpretation techniques (see Appendix G), the 
recommended curve for clean soils is characterized with an average probability of 26 %. 


Caution should be exercised when applying the procedure to sites where conditions are 
different from the case history data. The case history data used to develop the procedure are 
limited to relatively level ground sites with the following general characteristics: (1) uncemented 
soils of Holocene age; (2) average depths less than about 10 m; and (3) ground water table 
depths between 0.5 m and 6m. All Vs measurements are from below the water table. About 
three-quarters of the case history data are for soils with fines content greater than 5 %. Almost 
half of the case histories are for earthquakes with magnitudes near 7. 


Three concerns when using shear wave velocity as an indicator of liquefaction resistance 
are (1) its higher sensitivity (when compared with the penetration-based methods) to weak 
interparticle bonding, (2) the lack of a physical sample for identifying non-liquefiable clayey 
soils, and (3) not detecting thin liquefiable strata because the test interval is too large. The 
preferred practice is to drill sufficient boreholes and conduct sufficient other in situ tests to 
detect thin liquefiable strata, identify non-liquefiable clay-rich soils, identify silty soils above the 
ground water table that might have lower values of Vs should the water table rise, and detect 
liquefiable weakly cemented soils. 


4.2 FUTURE STUDIES 
The following future studies are recommended: 


1. Additional well-documented case histories with all types of soil that have and have 
not liquefied during earthquakes should be compiled, particularly from deeper deposits (depth > 
8 m) and from denser soils (Vs > 200 m/s) shaken by stronger ground motions (Ama > 0.4 g), to 
further validate the recommended curves. Also, case histories from lower magnitude 
earthquakes (MM, < 7) may improve estimates of the magnitude scaling factor. 


2. Laboratory and field studies should be conducted to further refine estimates of /5,, 
the limiting value of Vs; for cyclic liquefaction occurrence. For example, careful laboratory 
studies may identify more clearly the influence of fines content, gravel content, and particle 
gradation on )’s,. Additional careful penetration-Vs correlation studies may also help refine the 
Vs, estimates. 
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3. Laboratory studies should also be conducted to evaluate the implied assumption 
observed in Fig. 2.3 that at low values of Vs, (say 100 m/s) liquefaction resistance is independent 
of fines content. 


4. Additional work is needed to evaluate the significance of ignoring soil type and 
horizontal stress in the overburden correction. 

5. Standard test procedures exist only for the crosshole test. Standard test methods 
should be developed for the other im situ seismic tests. 
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APPENDIX B 


SYMBOLS AND NOTATION 


The following symbols and notation are used in this report: 


A 

a 

a}, A2, a3 
Qmax 

B,, B2 

b 

b;, b2, bs, be 
CRR 


KEmin) 
KVs1) 
I) 

Kiyo) 
jus) 


Inu(F's) 


parameter that depends on soil structure; 

parameter related to slope of CRR-Vs; curve; 

regression coefficients; 

peak horizontal ground surface acceleration; 

parameters relating Vs; and penetration resistance; 

parameter related to slope of CRR-Vs, curve; 

regression coefficients; 

average cyclic resistance ratio; 

CRR for cyclic triaxial tests; 

CRR corrected for high overburden stress; 

CRR for magnitude 7.5 earthquakes; 

cyclic stress ratio; 

factor to correct Vs for overburden pressure; 

median grain size by mass; 

double-amplitude axial strain; 

relative density; 

the constant e raised to the power of a given number; 

age and soil type factors for correlating Vs and N,; 

fines content (particles smaller than 75 pm); 

factor of safety; 

high overburden stress exponent; 

function of minimum void ratio; 

function of Vs1; 

function of x (= Fs); 

function of average peak cyclic shear strain; 

probability density function of calculated F's for 
liquefaction case histories; 

probability density function of calculated F's for 
non-liquefaction case histories; 
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shear modulus; 

small-strain shear modulus; 

Gmax corrected for confining stress and void ratio; 
secant shear modulus at Yay; 

acceleration of gravity; 

cementation and aging correction factor; 

fines content correction factor; 

coefficient of effective lateral earth pressure at rest; 
high overburden stress correction factor; 

natural logarithm function; 

magnitude scaling factor; 

earthquake moment magnitude; 

stress exponent; 

SPT blow count in Japanese practice; 

SPT energy-corrected blow count; 

SPT energy- and overburden stress-corrected blow count; 
magnitude scaling factor exponent; 

reference overburden stress (= 100 kPa); 

probability of liquefaction occurrence 

factor to account for effects of multidirectional shaking; 
shear stress reduction coefficient; 

sine function; 

residual standard deviation; 

small-strain shear wave velocity; 

overburden stress-corrected Vs; 

equivalent clean soil value of Vs;; 

limiting upper value of Vs; for liquefaction occurrence; 
mean stress-corrected Vs; 

depth; 

function of depth; 

function of depth; 

average peak cyclic shear strain; 

mean; 

mass density of soil; 

standard deviation; 

cyclic deviator stress in cyclic triaxial tests; 

initial effective horizontal confining stress; 

mean effective confining stress; 

initial effective confining stress in cyclic triaxial tests; 
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Oy 


Oy 


Tav 


Tmax 


total vertical (or overburden) stress; 

initial effective vertical (or overburden) stress; 

average cyclic equivalent uniform shear stress generated by 
earthquake; and 

maximum cyclic shear stress generated by earthquake. 
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APPENDIX C 


GLOSSARY OF TERMS 


The following definitions apply to this report: 


Case History 


Cnitical Layer 


Liquefaction 
Occurrence 

Moment Magnitude 
Overburden Stress- 
Corrected Shear Wave 
Velocity 

Peak Horizontal Ground 


Surface Acceleration 


Shear Wave Velocity 


Shear Wave 


An earthquake and a test array. 


The layer of non-plastic soil below the ground water table where 
corrected values of shear wave velocity and penetration are the 
least, and where cyclic stress ratios are the greatest. 


Surface manifestations of excess pore-water pressure at depth, 
such as sand boils, ground cracks and fissures, and ground 
settlement. 


An earthquake magnitude scale defined in terms of energy. 
Shear wave velocity measurement corrected to a reference 


vertical (or overburden) stress of 100 kPa. 


The peak value in a horizontal ground surface acceleration 
record that would occur at the site in the absence of liquefaction 
Or excess pore-water pressures. 


The velocity of a propagating shear wave within a material with 
either the direction of wave propagation or the direction of 


particle motion in the vertical direction. 


A body wave with the direction of particle motion transverse to 


~ the direction of wave propagation. 
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Test Array 


The two boreholes used for crosshole measurements, the borehole 
and source used for downhole measurements, the cone sounding 
and source used for seismic cone measurements, the borehole 
used for suspension logging measurements, or the line of 
receivers used for Spectral-Analysis-of-Surface-Waves (SASW) 
measurements. 
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APPENDIX D 


COMPARISON OF Vs-BASED LIQUEFACTION RESISTANCE CURVES 


During the past two decades, several studies have been conducted to investigate the 
relationship between Vs and liquefaction resistance. These studies involved laboratory tests 
(Dobry et al., 1981; Dobry et al., 1982; de Alba et al., 1984; Hynes, 1988; Tokimatsu and 
Uchida, 1990; Tokimatsu et al., 1991a; Rashidian, 1995; Rauch et al., 2000), analytical 
investigations (Bierschwale and Stokoe, 1984; Stokoe et al., 1988c; Andrus, 1994), penetration- 
Vs correlations (Seed et al., 1983; Lodge, 1994; Kayabali, 1996; Rollins et al., 1998b; Andrus et 
al., 1999), or field performance observations (Stokoe and Nazarian, 1985; Robertson et al., 
1992; Kayen et al., 1992; Andrus and Stokoe, 1997; Andmus et al., 1999; Juang and Chen, 2000; 
Andrus and Stokoe, 2000; Juang et al., 2001a). Many of the liquefaction evaluation procedures 
developed from these studies follow the general format of the Seed-Idriss simplified procedure, 
where Vs is corrected to a reference overburden stress and correlated with the cyclic stress, or 
resistance, ratio. 


This appendix reviews seven proposed liquefaction evaluation curves based on CRR and 
Vs;. The seven CRR-Vs; curves are shown in Fig. D.1. Each of the curves is briefly discussed 
below. 


D.1 CURVE BY TOKIMATSU AND UCHIDA (1990) 


The best-fit curve by Tokimatsu and Uchida (1990) shown in Fig. D.1 was determined 
from laboratory cyclic triaxial test results for various sands with less than 10 % fines (silt and 
clay) and 15 cycles of loading. Figure D.2 presents the cyclic triaxial test results. The solid 
symbols in Fig. D.2 correspond to specimens obtained by the in situ freezing technique. The 
open symbols correspond to specimens reconstituted in the laboratory. Tokimatsu and Uchida 
defined the cyclic resistance ratio for cyclic triaxial tests, CRR., as the ratio of cyclic deviator 
stress to initial effective confining stress, ¢g/2o', when the double-amplitude (or peak-to-peak) 
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*Curve adjusted using scaling 
factor of 1.19 for magnitude 
7 earthquakes 
*“Approximate curve for clean 
sand & 15 cycles of loading, 
assuming @mjn = 0.65, 
K, = 0.5, r- = 0.9 


Tokimatsu & 
Uchida (1990)** 


Liquefaction 


Robertson et al. 
(1992) 


Cyclic Resistance Ratio, CRR 


Liquefaction 


0 100 200 300 


Overburden Stress-Corrected Shear Wave 
Velocity, Vs7, m/s 


Fig. D.1 - Comparison of Seven Proposed CRR-Vs; Curves for Clean Granular Soils. 
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Fig. D.2 - Relationship Between Liquefaction Resistance and Normalized Shear Modulus 
for Various Sands with Less than 10 % Fines Determined by Cyclic Triaxial 
Testing. (modified from Tokimatsu and Uchida, 1990) 
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axial strain, DA, reaches 5 %. They measured the elastic shear modulus of the specimen at a 
shear strain of 10° % just prior to the liquefaction test. This small-strain shear modulus was 
normalized to correct for the influence of confining pressure and void ratio by: 


Gmax 
CS (D.1) 
Sem owe. 
and 
ag Ik 
flemin) = otal (D.2) 
+ @€min 
where 
Gy | = _ the normalized shear modulus, 
€mn = the minimum void ratio determined by standard test method, and 


O'm the mean effective confining stress. 

Tokimatsu and Uchida selected an exponent of 2/3 rather than 1/2, as determined by Hardin and 
Dmevich (1972), because it seemed that a slightly better correlation could be obtained. Values 
of €min ranged from 0.61 to 0.91 for the sands tested. The actual values of void ratio in each test 
were greater than @min, With values ranging from about 0.65 to about 1.4. 


By combining Egs. (1.1) and (D.1), one obtains the following relationship for converting 
Gy to mean stress-corrected Vs: 


0.33 0.5 
1 
Vsim = re | = [Gx Lean) (D.3) 
os fe) 
where 
Vsim = Mean stress-corrected Vs, and 
o’, = the mean effective confining stress in kgf/cm? (1 kgf/cm? = 98 kPa). 


Tokimatsu and Uchida (1990) suggested using 0.65 as an average value Of @min for clean sands. 
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The overburden stress-corrected Vs and Vs; can be related by: 


0.33 0.33 0.08 0.33 
] 3 ] 3 
V7 m — V: ~ V. aie al 4 
4 {4 te) ‘5 (= Gia Ve 


where 
Ke = ~the coefficient of lateral earth pressure at rest (=o';,/o', ). 


Values of Vs; for the best fit curve by Tokimatsu and Uchida (1990) shown in Fig. D.1 are 
determined (Andrus et al., 1999; after Tokimatsu et al., 1991a) from Fig. D.2 using Eqs. (D.3) 
and (D.4), and assuming K’, = 0.5, min = 0.65, o', = 100 kPa, and soil density of 1.9 Mg/m’. 


For converting CRR, to an equivalent field CRR, Tokimatsu and Uchida (1990) 
suggested the following expression originally proposed by Seed (1979): 


CRR = foeke 


r- (CRRx) | (D.5) 


where 


fe = a constant to account for the effects of multi-directional shaking with a value 
between 0.9 and 1.0. 


Values of CRR for the best fit curve by Tokimatsu and Uchida shown in Fig. D.1 are determined 
from Fig. D.2 using Eq. (D.5) and assuming K’, = 0.5 and r, = 0.9. 


Because the other liquefaction resistance curves shown in Fig. D.1 were drawn to bound 
liquefaction case histories, the more conservative lower bound curve for the laboratory test 
results by Tokimatsu and Uchida (1990) also is shown. This curve was drawn (Andrus et al., 
1999) from Fig. D.2 following the procedure outlined above. 


D.2 CURVE BY ROBERTSON ET AL. (1992) 


The bounding curve by Robertson et al. (1992) was developed using field performance 
data from primarily sites in Imperial Valley, California, along with data from four other sites, as 
shown in Fig. D.3. The soil at these sites contained as much as 35 % fines. Robertson et al. 
corrected Vs using Eq. (2.7). The shape of their curve was based on the analytical results of 
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Fig. D.3 - Liquefaction Resistance Curve for Magnitude 7.5 Earthquakes and Case 
History Data from Robertson et al. (1992). 
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Bierschwale and Stokoe (1984). They reasoned that the curve should pass close to the Imperial 
Valley (Wildlife site) data point, since liquefaction did and did not occur at this site during the 
1987 Superstition Hills (44, = 6.5) and Elmore Ranch (M, = 6.2) earthquakes, respectively. 
Robertson et al. used magnitude scaling factors similar to those suggested by Seed and Idriss 
(1982), Column 2 of Table 2.2, to position their curve for magnitude 7.5 earthquakes. 


D.3 CURVE BY KAYEN ET AL. (1992) 


Kayen et al. (1992) studied four sites that did and did not liquefy during the 1989 Loma 
Prieta, California, earthquake (M4, = 7.0). The four sites are: Port of Richmond, Bay Bridge 
Toll Plaza, Port of Oakland, and Alameda Bay Farm Island South Loop Road. The fines 
content for soils at these sites ranged from less than 5 % to as much as 57 %. Values of Vs were 
measured by the SCPT method and corrected for overburden stress using Eq. (2.7). Figure D.4 
presents their data and bounding curve. The curve by Kayen et al. shown in Fig. D.1 was 


adjusted for magnitude 7.5 earthquakes by assuming a MSF of 1.19 (see Column 3 of Table 
22): 


D.4 CURVE BY LODGE (1994) 


Lodge (1994) considered the same sites that Kayen et al. (1992) studied, as well as other 
sites shaken by the 1989 Loma Pneta earthquake. The curve by Lodge was developed as 
follows. First, cyclic stress ratios for the entire soil profile at each site were calculated. Second, 
available SPT blow counts were corrected for overburden pressure and energy. Soil layers with 
high and low liquefaction potential were identified with the procedure of Seed et al. (1985). Soil 
layers with corrected blow count within 3 of the SPT-based curve were eliminated due to 
uncertainties in the correlation. Third, V; measurements from SCPT and crosshole tests were 
corrected for overburden stress using Eq. (2.7). Fourth, on a “meter by meter” basis, values of 
Vs; and cyclic stress ratio were plotted for both layer types, those which were predicted 
liquefiable and those which were predicted non-liquefiable. Fifth, published data for sites shaken 
by the 1983 Borah Peak, Idaho, and 1964 Niigata, Japan, earthquakes were added to the plot. 
Finally, a curve was drawn to include all liquefiable layers, as shown in Fig. D.5. The curve by 
Lodge shown in Fig. D.1 was adjusted for magnitude 7.5 earthquakes by assuming a MSF of 
1.19 (see Column 3 of Table 2.2). 
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Fig. D.4 - Liquefaction Resistance Curve for Magnitude 7 Earthquake and Case History 
Data from Kayen et al. (1992). 
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Fig. D.5 - Liquefaction Resistance Curve for Magnitude 7 Earthquakes and Case 
History Data from Lodge (1994). 
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D.5 CURVE BY ANDRUS AND STOKOE (1997) 


The curve by Andrus and Stokoe (1997) shown in Fig. D.1 was developed for the 
proceedings of the 1996 NCEER Workshop (Youd and Idriss, eds., 1997). Several suggestions 
were offered at, and after, the workshop concerning how site variables should be define, as well 
as the shape of the boundary curve separating liquefaction and no liquefaction. Following the 
suggestions and using field performance data from 20 earthquakes and in situ Vs measurements 
from over 50 sites in soils ranging from clean fine sand to sandy gravel with cobbles to profiles 
including silty clay layers, Andrus and Stokoe constructed curves for uncemented, Holocene-age 
soils with various fines content. The values of Vs were corrected using Eq. (2.7). The curve by 
Andrus and Stokoe (1997) for fines content < 5 % along with the case history data are presented 
in Fig. D.6. The shape of the curve by Andrus and Stokoe (1997) was based on a modified 
relationship between Vs; and CSR for constant average cyclic shear strain suggested by R. 
Dobry. 


D.5.1 Cyclic Shear Strain and Cyclic Shear Stress 


Liquefaction results from the rearranging of soil particles and the tendency for decrease 
in volume. Experimental and theoretical studies show that decrease in volume is more closely 
related to cyclic strain than cyclic stress (Silver and Seed, 1971); a threshold cyclic strain exists 
below which neither rearrangement of soil particles nor decrease in volume take place (Dmevich 
and Richart, 1970; Youd, 1972; Pyke et al., 1975), and no pore water pressure buildup occurs 
(Dobry et al., 1981; Seed et al., 1983); and that there is a predictable correlation between cyclic 
shear strain and pore pressure buildup of saturated soils (Martin et al., 1975; Park and Silver, 
1975; Finn and Bhatia, 1981; Dobry et al., 1982; Hynes, 1988). The threshold cyclic strain is 
limited to a narrow range of variation, ranging from about 0.005 % for gravels to 0.01 % for 
normally consolidated clean sands and silty sands to 0.03 % for overconsolidated clean sands. 
In addition, cyclic strain-controlled test results are less affected than stress-controlled tests by 
factors such as density, confining stress, anisotropic confining stress, fabric and prestaining 
(Martin et al., 1975; Dobry and Ladd, 1980; Dobry et al., 1982; Hynes, 1988). It should also be 
noted that the steady state approach to liquefaction evaluation by Poulos et al. (1985) is based 
on a triggering strain level. These findings confirm the fact that cyclic strain is more 
fundamentally related to pore pressure buildup than cyclic stress, and are strong arguments in 
favor of a cyclic strain approach to liquefaction evaluation. 
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Fig. D.6 - Liquefaction Resistance Curve for Magnitude 7.5 Earthquakes and 
Uncemented Clean Soils of Holocene Age with Case History Data from 
Andrus and Stokoe (1997). 


Ie: 


Cyclic shear strain and cyclic shear stress can be related by the following equation: 


(D.6) 


where 


Vow = the average peak cyclic shear strain during a cyclic stress-controlled test of 
_ uniform cyclic shear stress Za, and 
(G) = the secant shear modulus at 7,, during the same cyclic test. 


In the cyclic strain approach proposed by Dobry et al. (1982), the average cyclic shear 
strain caused by an earthquake is estimated from: 


Amax Ovl'd Gmax 
ey SONG See te eel Come (D.7) 
& pVvs G),_ 


Equation (D.7) is obtained by combining Eqs. (1.1), (2.1) and (D.6). The vanation of shear 
modulus with strain is commonly expressed in terms of (G),_/Gmax called the modulus 
reduction factor. The modulus reduction factor can be estimated from an experimentally 
determined correlation. Neither pore pressure buildup nor liquefaction will occur when 7g, is 
less than the threshold strain. When 7,, is greater than the threshold strain, then pore pressure 
buildup can occur. The amount of pore pressure buildup can also be estimated from an 
experimentally determined correlation. 


D.5.2 Dobry’s Relationship Between CRR and Vs; 
R. Dobry (personal communication to R. D. Andrus, 1996) derived a relationship 


between Vs; and CSR for constant average cyclic shear strain using Eqs. (1.1) and (D.6). 
Combining Eqs. (1.1) and (D.6), and dividing both sides by o’, leads to: 


G 
Tav aro 2) G),. V3 (D.8) 
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For an overburden stress of 100 kPa, Vs = Vs, and curves of constant average cyclic strain can 
be expressed by: 


CSR = ms = fete) Vs)? (D.9) 
where 
Kya) = to £) ),, (D.10) 


Since CSR equals CRR at the point separating liquefaction from no liquefaction, Eq. (D.9) 
provides an analytical basis for establishing the CRR-Vs; curve at low values of Vs; (say Vs; < 
125 m/s) and extending them to zero at Vs; = 0. 


D.5.3 Modified CRR-Vs; Relationship 


Andrus and Stokoe (1997) reasoned that the curve separating liquefiable and non- 
liquefiable soils would become asymptotic to some limiting upper value of Vs;._ They modified 


Eq. (D.9) to: 


Vs) 1 i 
CRR = of 22) if : MSF 11 
| 100 Vs Visie Visi ay 
where 
Vs, = the limiting upper value of Vs, for liquefaction occurrence, and 
a,b = curve fitting parameters. 


The first term in Eq. (D.11) is a form of Eq. (D.9), assuming (yay) is independent of initial 
effective confining pressure and of pore water pressure buildup. The second term is a hyperbola 
with a small value at low values of Vs;, and a very large value as Vs; approaches )’5,. 


The curve by Andrus and Stokoe (1997) shown in Figs. D.1 and D.6 is defined by Eqs. 
(D.11) and (2.9) with a = 0.03, b= 0.9, n = -3.3, andV’s, = 220 ms. 


1 


D.6 CURVE BY ANDRUS ET AL. (1999) 


Since the publication of the 1996 NCEER Workshop proceedings (Youd and Idriss, 
eds., 1997), the case history data compiled by Andrus and Stokoe (1997) have been revised 
based on new information, and expanded to include field performance data from 26 earthquakes 
and more than 70 measurements sites. Also, the 1998 MCEER Workshop was held to discuss 
developments since the 1996 workshop. From the suggestions given at the second workshop 
and using the expanded database, the curve proposed by Andrus and Stokoe (1997) was revised 
in the report by Andrus et al. (1999) and paper by Andrus and Stokoe (2000). The revised 
curve for uncemented soils with fines content < 5 % along with the case history data are shown 
in Fig. D.7. The development of the revised curve is discussed in Appendix F. 


D.7 SUMMARY 


Seven proposed curves relating CRR and Vs; were discussed in this Appendix. Many of 
the differences among the seven curves (see Fig. D.1) can be explained by the following four 
factors: (1) The best-fit curve by Tokimatsu and Uchida (1990) is more of a median curve, 
while the other curves bound the liquefaction case history data. (2) Portions of the proposed 
curves are based on limited data, and the investigator(s) have assumed different levels of 
conservatism. In particular, the curves by Robertson et al. (1992), Kayen et al. (1992), and 
Lodge (1994) were based on little or no data above Vs; of 200 m/s, and were conservatively 
drawn in this region. (3) Methods for selecting some site variables and correction factors are 
different among investigator(s). (4) Some errors exist in the database by Andrus and Stokoe 
(1997), and lead to more conservative curve than the updated curve by Andrus et al. (1999) 
above a Vs; value of 150 m/s. Thus, the CRR-Vs; curve proposed by Andrus et al. (1999) for 
clean soils is recommended because it was based on the largest, most correct case history data 
set and procedures recommended by the 1996 NCEER Workshop (Youd et al., 1997). 
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Fig. D.7 - Revised Liquefaction Resistance Curve for Magnitude 7.5 Earthquakes and 
Uncemented Clean Soils of Holocene Age with Case History Data from 
Andrus et al. (1999). 
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CASE HISTORY DATA AND THEIR CHARACTERISTICS 
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Shear wave velocity measurements have been made for field liquefaction studies at many 
sites during the past twenty years. Table E.1 presents a list of over 70 sites and 26 earthquakes 
that have been investigated. Of the 26 earthquakes listed, 9 occurred in the United States; and 
the other 15 in Japan, Taiwan, and China. The field performance information for these 
earthquakes along with the Vs measurements provides an important opportunity to determine the 
relationship between liquefaction resistance and Vs directly from case histories. A detailed 
summary of available case history data is presented in Appendix H. This appendix describes the 
site variables and characteristics of the database. 


E.1 SITE VARIABLES AND DATABASE CHARACTERISTICS 
E.1.1 Earthquake Magnitude 


Earthquake magnitudes for the 26 earthquakes listed in Table E.1 range from 5.3 to 8.3, 
based on the moment magnitude scale. Moment magnitude is the scale most commonly used for 
engineering applications, and is the preferred scale for liquefaction resistance calculations (Youd 
et al., 1997). When other magnitude scales are reported by the investigator(s), they are 
converted to M,, using the relationship of Heaton et al. (1982) shown in Fig. 2.5. 


E.1.2 Shear Wave Velocity Measurement 


At the more than 70 investigation sites listed in Table E.1, shear wave velocity 
measurements were made with 139 test arrays. A test array is defined in this report as the two 
boreholes used for crosshole measurements, the borehole and source used for downhole 
measurements, the cone sounding and source used for seismic cone measurements, the borehole 
used for suspension logging measurements, or the line of receivers used for Spectral-Analysis- 
of-Surface-Waves (SASW) measurements. Of the 139 test arrays, 39 are crosshole, 21 
downhole, 27 seismic cone, 15 suspension logger, 36 SASW, and one is unknown. 
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Table E.1 - Earthquakes and Sites Used to Establish Liquefaction Resistance Curves 


Earthquake Moment Site Reference 
Magnitude 
(1) (2) (3) (4) 
1906 San Francisco, California es Coyote Creek; Salinas River | Youd & Hoose (1978); 
(North, South) Barrow (1983); Bennett & 
Tinsley (1995) 


Marina District (2, 3, 4, 5, 
School) 


1957 Daly City, California of) Kayen et al. (1990); 
Tokimatsu et al. (1991b); 
T. L. Youd (personal 
communication to R. D. 
Andrus, 1999) 
1964 Niigata, Japan 
ie 


1975 Haicheng, China 3 


Niigata City (Al, Cl, C2, 
Railway Station) 


Yoshimi et al. (1984; 
1989); Tokimatsu et al. 
(1991a) 


Chemical Fiber; 
Construction Building; 
Fishery & Shipbuilding; 
Glass Fiber; Middle School; 


Arulanandan et al. (1986) 


1979 Imperial Valley, California 
1981 Westmorland, California 
1987 Elmore Ranch, California 
1987 Superstition Hills, California 


Heber Road (Channel fill, 
Point bar); Kornbloom; 
McKim; Radio Tower; Vail 
Canal; Wildlife 


Bennett et al. (1981; 
1984); Sykora & Stokoe 
(1982); Youd & Bennett 
(1983); Bierschwale & 
Stokoe (1984); Stokoe & 
Nazarian (1984); Dobry 
et al. (1992); Youd & 
Holzer (1994) 


1980 Mid-Chiba, Japan 5.9 Owi Island No. 1 Ishihara et al. (1981; 
1985 Chiba-Ibaragi-Kenkyo, Japan 6.0 1987) 


1983 Borah Peak, Idaho Andersen Bar; Goddard Youd et al. (1985); Stokoe 
Ranch; Mackay Dam et al. (1988a); Andrus et 

Downstream Toe; North al. (1992); Andrus (1994) 
Gravel Bar; Pence Ranch 


1986 Event LSST2, Taiwan 
Event LSST3, Taiwan 
Event LSST4, Taiwan 
Event LSST6, Taiwan 
Event LSST7, Taiwan 
Event LSST8, Taiwan 
Event LSST12, Taiwan 
Event LSST13, Taiwan 
Event LSST16, Taiwan 


1987 Chiba-Toho-Oki, Japan | 65 | Sunamachi Ishihara et al. (1989) 
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Lotung LSST Facility Shen et al. (1991); 


EPRI (1992) 


Table E.1 (cont.) - Earthquakes and Sites Used to Establish Liquefaction Resistance Curves. 


Earthquake 


(1) 


1989 Loma Prieta, California 


1993 Kushiro-Oki, Japan 


1993 Hokkaido-Nansei-Oki, Japan 


1994 Northridge, California 


Site 


Bay Bridge Toll Plaza, Bay 
Farm Island (Dike, South 
Loop Road); Port of 
Oakland; Port of Richmond 


Coyote Creek; Salinas River 
(North, South); 


Marina District (2, 3, 4, 5, 
school) 


Moss Landing (Harbor 
Office, Sandholdt Road, 
State Beach) 


Santa Cruz (SC02, SC03, 
SC04, SC05, SC13, SC14) 


Treasure Island Fire Station 


Treasure Island Perimeter 
(Approach to Pier, UMO3, 
UM05, UM06, UM09) 


Kushiro Port (2, D) 


Pension House 


Hakodate Port 
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Reference 


(4) 


Stokoe et al. (1992); 
Mitchell et al. (1994) 


Barrow (1983); 

M. J. Bennett (personal 
communication to R. D. 
Andrus, 1995); Bennett 
and Tinsley (1995) 


Kayen et al. (1990); 
Tokimatsu et al. (1991b) 


Boulanger et al. (1995); 
Boulanger et al. (1997) 


Hryciw (1991); 
Hryciw et al. (1998) 


Hryciw et al. (1991); 
Redpath (1991); Gibbs 
et al. (1992); Furhriman 
(1993); Andrus (1994); 
de Alba et al. (1994) 


Geomatrix Consultants 
(1990); Hryciw (1991); 

R. D. Hryciw (personal 
communication to R. D. 
Andrus, 1998); Hryciw 

et al. (1998); Andrus et al. 
(1998a, 1998b) 


Jai et al. (1995); S. Iai 
(personal communication 
to R. D. Andrus, 1997) 


Kokusho et al. (1995a, 
1995b, 1995c) 


S. Iai (personal 
communication to 
R. D. Andms, 1997) 


Abdel-Haq & Hryciw 
(1998) 


Table E.1 (cont.) - Earthquakes and Sites Used to Establish Liquefaction Resistance Curves. 


Earthquake Moment Site Reference 
Magnitude 
(1) (2) (3) (4) 


1995 Hyogo-Ken Nanbu, Japan Hanshin Expressway 5 
(3, 10, 14, 25, 29); Kobe- 
Nishinomiya Expressway 
(3, 17, 23, 28) 


KNK;; Port Island (Downhole 
Array); SGK 


Port Island (Common 
Factory) 


Kobe Port (7C); Port Island 
(1C, 2C) 


Kobe Port (LPG Tank Yard) 


Hamada et al. (1995); 
Hanshin Expressway 
Public Corporation (1998) 


Sato et al. (1996); 
Shibata et al. (1996) 


Ishihara et al. (1997); 
Ishihara et al. (1998) 


Inatomi et al. (1997); 
Hamada et al. (1995) 


S. Yasuda (personal 
communication to 
R. D. Andrus, 1997) 


Values of Vs reported by the investigator(s) are used directly. The one exception is for 
the downhole array located at the Marina District School site in San Francisco, California. A 
reevaluation of the field data indicates that Vs values reported for the critical layer at this site are 
too high. They are recalculated using the pseudo-interval method, as discussed in Section E.2.2. 
Only the crosshole measurements made with shear waves having particle motion in the vertical 
direction are used. Crosshole measurements near the critical layer boundary that seem high, and 
could represent refracted waves, are not included in the average. Some Vs values are from 
measurements made before the earthquake, others followed the earthquake. No adjustments are 
made to compensate for changes in soil density and Vs due to ground shaking. 


E.1.3 Measurement Depth. 


In situ Vs measurements may be reported at discrete depths or for continuous intervals, 
depending on the test method. When velocities are reported for continuous intervals, as is 
typically the case for downhole, seismic cone, suspension logger and SASW measurements, the 
depth to the center of each interval is assumed. Thus, if the reported Vs profile has ten velocity 


layers, it is assumed that the profile consists of ten “measurements” with depths at the center of 
each layer. 
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E.1.4 Case History 


In this report, a case history is defined as a seismic event and a test array. For example, 
at the Treasure Island Fire Station site, crosshole measurements were made between five 
different pairs of boreholes, downhole measurements were made by two different investigators, 
seismic cone measurements were made at one location, and SASW measurements were made 
along one alignment. Thus, a total of nine case histories are identified for the Fire Station site 
and the 1989 Loma Prieta, California earthquake. At the Marina District School site, downhole 
measurements were made at one location. Estimates of ground surface acceleration at this site 
are available for the 1957 Daly City and 1989 Loma Pneta earthquakes. Thus, two case 
histories are identified for the Marina District School site. Combining the 26 seismic events and 
139 test arrays, a total of 225 case histories are obtained with 149 from the United States, 36 
from Taiwan, 34 from Japan, and 6 from China. 


The two exceptions to this definition are the Owi Island No. 1 site and the Moss Landing 
Sandholdt Road UC-4 site where additional subsurface information is available. At Owi Island, 
pore pressure transducers recorded pore-water pressure buildup for two separate layers. At 
Moss Landing, inclinometer measurements indicated lateral movement in an upper loose layer 
and no lateral movement in a lower dense layer. Thus, two case histories are identified for each 
of these two test arrays. 


E.1.5 Liquefaction Occurrence 


It is important to realize that the occurrence of liquefaction, in this evaluation, is based 
on the appearance of surface evidence, such as sand boils, ground cracks and fissures, and 
ground settlement. Case histories are classified as non-liquefaction when no liquefaction effects 
were observed. At the Owi Island No. 1, Lotung LSST Facility, Sunamachi, Wildlife (1987 
earthquakes), and Port Island sites, the assessment of liquefaction or non-liquefaction 
occurrence is supported by pore-water pressure measurements. In addition, liquefaction 
occurrence is assigned (in This Report) to the Treasure Island, California, Fire Station case 
histories where the strong ground motion records from the 1989 Loma Prieta earthquake exhibit 
a sudden drop at about 15 seconds and small motion afterward (Idriss, 1990), indicating 
liquefaction (de Alba et al., 1994). Of the 225 case histories, 99 are liquefaction case histories 
and 126 are non-liquefaction case histories. Figure E.1 shows the distribution of case histories 
with earthquake magnitude. 
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Fig. E.1 - Distribution of Liquefaction and Non-Liquefaction Case Histories by Earthquake 
Magnitude. 


E.1.6 Critical Layer 


The layer of soil most likely to liquefy at a site, or the critical layer, is the layer of non- 
plastic soil below the ground water table where values of Vs;, as Gefined in Chavis 2, and 
penetration resistance are generally the least and cyclic stress Fauo relative to Vs; is the greatest. 
Figure E.2 presents the cumulative relative frequency distributions for we case histones by 
critical layer thickness and predominate soil type (gravel, or sand and silt). Critical layer 
thicknesses range from 1 mto as much as 15 m. About 50 % of the case niStGnes have a critical 
layer thickness less than 3.5 m; 90 % of the case histories have a critical layer thickness less than 
7m. Overall, the layer thicknesses for the gravel cases are less than the layer thicknesses for the 
sand and silt cases 


Figure E.3 presents the cumulative relative frequency distributions for the case histories 
by average Vs measurement depth in the critical layer and predominate oh type. The average 
measurement depths are between 2 m and 11 m for nearly all case histories. Over 50 % of the 
case histories have average measurement depths less than 5.5 m. About 90 % of the case 
histories have average measurement depths less than 8 m. Overall, the measurement depths for 
the gravel cases are less than the measurement depths for the sand and silt cases. 
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Fig. E.2 - Cumulative Relative Frequency of Case History Data by Critical Layer Thickness. 
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Fig. E.3 - Cumulative Relative Frequency of Case History Data by Average Depth of Vs 
Measurements in Critical Layer. 
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Materials comprising the critical layers range from clean fine sand to sandy gravel with 
cobbles to profiles including silty clay layers. In Fig. E.4, the distribution of case histories with 
earthquake magnitude, predominate soil type (gravel, or sand and silt) and average fines content 
(silt and clay) is presented. Of the 225 case histories, 28 were for sands with fines content (FC) 
< 5 %, 90 for sands with FC = 6 % to 34 %, 71 for sands and silts with FC > 35 %, 26 for 
gravels with FC < 5 %, and 10 for gravels with FC = 6 % to 34 %. 


About 70 % of the case histories are for natural soils deposits, with many formed by 
alluvial processes. The other 30 % are for hydraulic or dumped fills. Eight of the fills have been 
densified by soil improvement techniques. 


At least 85 % of the case histories are of Holocene age (< 10 000 years). Although the 
age of the other 15 % is unknown, they are believed to be also of Holocene age. 
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Fig. E.4 - Distribution of Case Histories by Earthquake Magnitude, Predominate Soil Type, and 
Average Fines Content. | 
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E.1.7 Ground Water Table 


Figure E.5 presents the cumulative relative frequency distnbutions for the case histories 
by depth to the ground water table and predominate soil type. The ground water table for nearly 
all case histories lies between depths of 0.5 m and 6 m. Nearly 60 % of the case histories have 
water table depths less than 2m. About 90 % of the case histories have water table depths less 
than 4.5 m. 


Artesian pressures are reported for the Lotung Large-Scale Seismic Test (LSST) Facility 
site in Tatwan. At this site, the pore-water pressure distribution is assumed to vary linearly from 
a pressure head of 8.1 m at a depth of 7 m to a pressure head of 1.9 m at a depth of 2 m. 


E.1.8 Total and Effective Overburden Stresses 
Values of total and effective overburden stresses are estimated using densities reported 
by the investigator(s). When no densities are reported, typical values for soils with similar grain 


size, penetration and velocity characteristics are assumed. In most instances, the assumed 
densities are 1.76 Mg/m? for soils above the water table and 1.92 Mg/m’ for soils below the 


water table. 
. 24 
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Fig. E.5 - Cumulative Relative Frequency of Case History Data by Depth to the Ground Water 
Table. 
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E.1.9 Average Peak Ground Acceleration 


Average values of peak horizontal ground surface acceleration, Qmax, are determined by 
averaging estimates reported by the investigator(s) and estimates made as part of this study 
using attenuation relationships developed from published ground surface acceleration data. 
Because many published attenuation relationships are based on both peak values obtained from 
ground motion records for the two horizontal directions (sometimes referred to as the randomly 
oriented horizontal component), the geometric mean (square root of the product) of the two 
peak values is used. Use of the geometric mean is consistent with the development of the SPT- 
based procedure (Youd et al., 1997; 2001). For the cases in this study, the difference between 
the geometric mean and arithmetic mean values is generally small, within about 5 %. 


E.1.10 Average Cyclic Stress Ratio 


Cyclic stress ratios, CSR, are first calculated for each “measurement” depth within the 
critical layer using Eq. 2.1 and then averaged. Values of rq are estimated using the average 
relationship developed by Seed and Idriss (1971) shown in Fig. 2.1. These 7z values are used to 
follow the traditional format of the SPT- and CPT-based procedures where the magnitude 
scaling factor is used to account for all effects of earthquake magnitude. 


E.1.11 Average Overburden Stress-Corrected Shear Wave Velocity 


Values of Vs within the critical layer are first corrected for overburden stress using Eq. 
2.7 and then averaged. The number of values included in the average range from 1 to 22 (see 
Appendix H). Values of Cys used to correct measured shear wave velocities range from 1.4 to 
0.9 for most of the data. About 80 % of the case histories have two to seven values in the 
average. No adjustments are made for possible variations between seismic test methods due to 
different source-receiver orientations with respect to the stress state in the soil. In the 
calculations, each site is assumed to be level ground. 
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E.2 SAMPLE CALCULATIONS 


Calculations for two sites shaken by the 1989 Loma Prieta, California earthquake (M,, = 
7.0) are presented below to illustrate how values of CSR and overburden stress-corrected shear 


wave velocity, Vs;, are determined. The two sites are Treasure Island Fire Station and Marina 
District School. 


E.2.1 Treasure Island Fire Station 


Treasure Island is a man-made island located in the San Francisco Bay along the Bay 
Bridge between the cities of San Francisco and Oakland. It was constructed in 1936-37 by 
hydraulic filling behind a perimeter rock dike. 


Extensive field tests have been conducted at the fire station on Treasure Island. Figure 
E.6 presents two Vs profiles for the site. The Vs profile determined by crosshole testing is from 
Fuhriman (1993). The other Vs profile is based on unpublished SASW test results by The 
University of Texas at Austin in 1992. Also presented in Fig. E.6 is the soil profile for the site. 
From the description by de Alba et al. (1994), the upper 4.5 m of soil consists of silty sand fill, 
possibly formed by dumping. Between depths of 4.5 m and 12.2 m, the soil consists of silty 
sand to clayey sand, formed by hydraulic filling. Beneath the hydraulic fill are natural clayey 
soils. The ground water table lies near the ground surface at a depth of 1.4 m. The critical layer 
is determined to be between depths of 4.5 m and 7 m, where the soil is non-plastic, lies below 
the water table, and exhibits the lowest values of Vs; relative to the highest values of CSR in the 
layer (see Fig. 3.1). 


During the 1989 Loma Prieta earthquake, a seismograph station at the fire station 
recorded ground surface accelerations. The peak values in the two horizontal accelerometer 
records are 0.16 g and 0.11 g (Brady and Shakal, 1994). 


Sample calculations for the crosshole and SASW test arrays are summarized in Tables 
E.2 and E.3, respectively. The data points used in the calculations are shown by the open 
symbols in Fig. E.6. Total and effective overburden stresses are calculated assuming densities of 
1.76 Mg/m? above the water table and 1.92 Mg/m’ below the water table. Stress reduction 
coefficients are estimated using the average curve by Seed and Idriss (1971) shown in Fig. 2.1. 
The geometric mean of the two peak values observed in the horizontal ground surface 
acceleration records is 0.13 g. Using these parameters, values of CSR and Vs; are calculated for 
the crosshole measurement at depth of 4.6 m as follows: 
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Fig. E.6 - Shear Wave Velocity and Soil Profiles for the Treasure Island Fire Station Site. 
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Representative values of CSR and Vs; used to defined the two case histories are determined by 
averaging values for each “measurement” depth within the critical layer, as shown in Tables E.2 
and E.3. 


E.2.2 Marina District School 


Kayen et al. (1990) conducted downhole seismic tests at the Winfield Scott School in the 
Marina District of San Francisco. Figure E.7 presents soil and velocity profiles for the site. The 
critical layer lies between depths of 2.7 m, the ground water table depth, and 4.3 m, the base of 
sand fill. The average Vs profile shown in Fig. E.7 was determined by Kayen et al., and was 
based on best-fit line segments through travel time measurements plotted versus depth. The 
second Vs profile is determined using the pseudo-interval method (This Report), as illustrated in 
Fig. E.8. Both methods should provide similar average values over the same depth interval. 
However, the layering assumed for the best-fit line segment method does not seem appropriate 
for the fill. For this reason, values of Vs based on the pseudo-interval method are used in this 
analysis. 


As discussed in Section 3.2, the Marina District of San Francisco experienced a peak 
horizontal ground surface acceleration of about 0.15 g during the 1989 Loma Pneta earthquake. 


Sample calculations for the Marina District School site are summarized in Table E.4. 
The locations of Vs measurements are assumed midway between receiver positions, as shown in 
Fig. E.7. Total and effective overburden stresses are estimated assuming densities of 
1.76 Mg/m? above the water table and 1.92 Mg/m’ below the water table. The ground water 
table is at a depth of about 2.7 m. Average values of CSR and Vs, defining the case history are 
determined by averaging values for the two “measurement” depths, as shown in Table E.4 
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Table E.2 - Sample Calculations for the Treasure Island Fire Station Site, Crosshole Test Array 
B1- B4, and the 1989 Loma Pneta Earthquake. 


Overburden 


Stress Stress- 
Effective Reduction Corrected 
Measurement | Average Overburden | Overburden | Coefficient’, Shear Wave 
Number Depth, Stress’, Stress), Velocity, Vs;, 
m kPa m/s 


(2) 


‘Assuming water table at 1.4 m; and material densities are 1.76 Mg/m’ above the water table 
and 1.92 Mg/m? below the water table. 

Based on average values determined by Seed and Idriss (1971). 

Assuming peak horizontal ground surface acceleration is 0.13 g. 


Table E.3 - Sample Calculations for the Treasure Island Fire Station Site, SASW Test Array, 
and the 1989 Loma Prieta Earthquake. 


Effective Reduction 


Measurement Overburden | Overburden | Coefficient’, Shear Wave 
Number Stress! , Stress!, rg Velocity, Vs;, 
m/s 


(6) 


and 1.92 Mg/m? below the water table. 
Based on average values determined by Seed and Idriss (1971). 
?Assuming peak horizontal ground surface acceleration is 0.13 g. 
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Fig. E.7 - Shear Wave Velocity and Soil Profiles for the Marina District School Site 
(Kayen et al., 1990). 
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Fig. E.8 - General Configuration of the Downhole Seismic Test Using the Pseudo-Interval 
Method to Calculate Shear Wave Velocity. 


E.3 SUMMARY 


The case history data described in this chapter are limited to level and gently sloping sites 
with the following characteristics: 


(1) average critical layer depths less than 10 m; 

(2) uncemented soils of Holocene age; 

(3) ground water table depths between 0.5 m and 6 m; and 
(4) all Vs measurements from below the water table. 


Of the 225 case histories, 57 are for soils with FC < 5 %, 98 for soils with FC = 6 % to 34 %, 


and 70 with FC > 35 %. About 20 % of the case histories are for soils containing more than 10 
% gravel. Nearly 50 % of the case histories are for earthquake magnitudes near 7. 
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Table E.4 - Sample Calculations for the Marina District School Site and the 1989 Loma Prieta 
Earthquake. 


Overburden 


Stress Stress- 
Total Effective Reduction Corrected 
Overburden | Overburden | Coefficient”, Shear Wave 
Stress’, Stress’, ra Velocity, 
kPa kPa Vs;, m/s 


(6) 


‘Based on pseudo-interval method. 

? Assuming water table at 2.7 m; and material densities are 1.76 Mg/m’ above the water table 
and 1.92 Mg/m? below the water table. 

3Based on average values determined by Seed and Idriss (1971). 

4Assuming peak horizontal ground surface acceleration is 0.15 g. 
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APPENDIX F 


DEVELOPMENT OF LIQUEFACTION RESISTANCE CURVES 
FROM CASE HISTORY DATA 


In the process of developing the liquefaction evaluation chart shown in Fig. 2.3 all case 
histories were initially plotted on the same chart. This aggregation was accomplished through 
an adjustment procedure; that is, the CSR values in each case history were adjusted to an 
earthquake with M, = 7.5 by dividing by Eq. (2.9) with m = -2.56. As done in penetration 
evaluation procedures, the sandy soil case histories were separated into three categories: (1) 
sands with average FC < 5 %; (2) sands with average FC = 6 % to 34 %; and (3) sands and silts 
with average FC > 35 %. For consistency, the gravelly soil case histones also were divided into 
the same three categories based on fines content. However, no case histories exist in the 
database with gravel having FC > 35 %. All data are plotted in Fig. F.1 along with the 
recommended CRR-Vs; curves. Development of these curves is discussed in this appendix. 


The shape of the CRR-Vs, curves shown in Fig. F.1 is based on a modified relationship 


between shear wave velocity and cyclic stress ratio for constant average cyclic shear strain 
suggested by R. Dobry. The modified relationship is expressed as (Andrus and Stokoe, 1997): 


2 
V si 1 1 
CRR75= o{ P| 14-5 | 1 
e 100 Psion yesie V'si ! 


where 


CRRz5= CRR for magnitude 7.5 earthquakes, 
ee the limiting upper value of Vs; for liquefaction occurrence, and 
a,b curve fitting parameters. 


As discussed in Section D.5, the first (quadratic) term of Eq. (F.1) is a form of Dobry’s 
relationship given by Eq. (D.9). The second term is a hyperbola with a small value at low values 
of Vs), and a very large value as Vs; approaches 5). 
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Fig. F.1 - Curves Recommended for Calculation of CRR from Shear Wave Velocity 
Measurements Along with Case History Data Based on Lower Bound 
Values of MSF for the Range Recommended by the 1996 NCEER 
Workshop (Youd et al., 1997) and rg Developed by Seed and Idriss (1971). 
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F.1 LIMITING UPPER Vs; VALUE FOR LIQUEFACTION OCCURRENCE 


As shown in Fig. F.1, CSR-values above about 0.35 are limited in the case history data. 
Thus, current estimates of )’5, rely, in part, on penetration-shear wave velocity correlations and, 
in part, on the data trend in Fig. F.1. 


F.1.1 Sandy Soils 


In the SPT-based procedure, a corrected blow count, (N;)60, of 30 is assumed as the 
limiting upper value for liquefaction occurrence in sands with < 5 % silt and clay (Seed et al., 
1985; Youd et al., 1997). Table F.1 presents estimates of equivalent Vs, for corrected blow 
count of 30. The correlation by Ohta and Goto (1978) modified to a blow count with a 
theoretical free-fall energy of 60 % (Seed et al., 1985) suggested equivalent Vs, values of 207 
m/s for Holocene sands, assuming that a depth of 10 m is equivalent to an effective overburden 
stress of 100 kPa. The stress-corrected crosshole measurements compiled by Sykora (1987b) 
for Holocene sands and non-plastic silty sands below the ground water table, with (Ni)60 
between 25 and 35, exhibited an average Vs; value of 206 m/s and standard deviation of 41 m/s. 
Finally, the case history data in this study were used to investigate the Vs; and (N1)60 relationship 
for well-documented sand layers with less than 10 % fines. These data are presented in Fig. F.2 
along with the best-fit relationship that can be expressed as: 


Vs; = B, [(v 1 Ie if (F.2) 


where B,; = 93.2 + 6.5 and B2 = 0.231 + 0.022 for soils with fines content < 10 %, and with Vs, 
in m/s and (Nj)60 in blows/0.3 m. The plotted data exhibit a mean Vs; value of 204 m/s at a 
(N,)6o Value of 30 and residual standard deviation, S,-;, of 12 m/s. 


In the CPT-based procedure, a normalized cone tip resistance, g.:v, of 160 is assumed as 
the limiting upper value for liquefaction occurrence in sands with < 5 % silt and clay (Youd et 
al., 1997; Robertson and Wnide, 1998). Figure F.3 presents average values of Vs; and q.in for 
soil layers with less than 10 % fines at several sites listed in Table E.1. Also shown in Fig. F.3 is 
the best-fit relationship for the plotted data, which can be expressed as: 


Vs, = B, [dew i (F.3) 
where B, = 88.2 + 15.5 and Bz = 0.154 + 0.037 for soils with fines content < 10 %, and with Vs; 
in m/s and q-:n is normalized tip resistance based on procedures by Robertson and Wride (1998). 


As noted in Table F.2, the plotted data exhibit a mean Vs; value of 193 m/s at a q.1n value of 160 
and residual standard deviation of 19 m/s. 
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Table F.1 - Estimates of Equivalent Vs, for Holocene Sands and Gravels Below the Ground 


Reference 
(1) 


Ohta & 
Goto 
(1978); also 
given in 
report by 
Sykora 
(1987a, 
page 29) 


Sykora 
(1987b, 
page 90); 
This Report 


Rollins et 
al. (1998a) 


This Report 
(see Fig. 
He) 


Water Table with Corrected SPT Blow Count of 30. 


Equivalent Vs) 
Estimate (m/s) 
(3) 


207 
... for Holocene 
sands 


Relationship 
(2) 
Vs = 69 (N) 2195 F,F> 
~ N; =SPT blow count 
measured in Japanese 


practice 

z =depth, m 227 

F, = 1.00 for Holocene- ...for Holocene 
age soils gravels 


F2 = 1.085 for sands; 
1.189 for gravel 
.. best-fit relationship for 289 
sets of SPT and Vs 
measurements from Ja 


Correlation between (N;)6o and 
crosshole Vs, normalized to 
effective overburden stress, 
measurements for Holocene 
sands and non-plastic silty 


206 
...for Holocene 
sands and non- 
plastic silty sands 
below the water 


sands below the ground water | table 
table at sites in U.S.A.; 16 sets | ...standard 
deviation is 


41 m/s 


V5= 53 (Neo)? (oy 232 


... best-fit relationship using ..for Neemahe 
equivalent Neo-values from ca bee f data Jj 
Becker Penetration Tests and ie 5 ee 


Vs measurements; 186 points 


from 7 Holocene gravel sites relationship 


he 204 
Vs) = Bi lw Deol? _..for Holocene 
B, =93.2+465 clean sands 


Bz = 0.231 + 0.022 
...best-fit relationship for 
uncemented, Holocene-age 


below the water 
table 


.. residual 
sands with less than 10 % non- standard 
plastic fines; 25 sets of average | qeviation is 
SPT and Vs; measurements all 12 m/s 


from below the water table 


100 


WN = 


WN 


Assumptions 
(4) 
. N; = 60/67 Neo 
‘ Neo = 30 


. Z= 10 mis equivalent to 


an overburden stress of 
100 kPa 


. All measurements are 


from below the ground 
water table 


. Average for Vs; values 
with (V;)6o between 


25 and 35 


. o'y = 100 kPa 


; Neo = 30 
. ov = 100 kPa 
. All measurements are 


from below the ground 
water table 


. Average for Vs; 


with (N1)6o = 30 


. ov = 100 kPa 
. Corrected blow count 


based on procedures 
given in Seed et al. 
(1985) and Robertson 
and Wride (1997; 1998) 
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Fig. F.2 - Variations in Vs, with (N;)69 for Uncemented, Holocene-age Sands with Less 
than 10 % Non-Plastic Fines. 
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Fig. F.3 - Variations in Vs, with g,;, for Uncemented, Holocene-age Sands with Less than 
10 % Non-Plastic Fines. 


102 


Table F.2 - Estimates of Equivalent Vs; for Holocene Sands Below the Ground Water Table 
with Normalized Cone Tip Resistance of 160. 


Equivalent Vs, 
Estimate (m/s) 
(3) 


193 
...for Holocene 
clean sands 
below the water 
table 


Reference 
(1) 


This Report 


(see Fig. 
F.3) 


Relationship 
(2) 


Vs; = B; Gan)” 
B, =88.2+15.5 
Bz =0.154 + 0.037 
..-best-fit relationship for 
uncemented, Holocene-age 


| 1. Average for Vs; 
with Wein = 160 

2. o'» = 100 kPa 

| 3. Normalized tip 

resistance based on 


sands with less than 10% non-__| ---"eSidual procedures given in 
plastic fines; 23 sets of average standard ; Robertson and Wride 
SPT and Vs measurements all deviation is (1997; 1998) 

19 m/s 


from below the water table 


From these estimates, a Vs; value of 210 m/s is assumed equivalent to an (Nz)eo value of 
30 in clean sands (< 5 % fines). A limiting upper Vs; value of 210 m/s for cyclic liquefaction 
occurrence at CSR = 0.6 is less than the general consensus value of 230 m/s suggested at the 
1998 MCEER Workshop. As a result, Figs. F.2 and F.3 were added specifically to provide 
additional evidence to support the use of 210 m/s in clean sands. 


For sandy and silty soils with FC > 35 %, the SPT-based chart by Seed et al. (1985) 
indicates a limiting upper (N)60 value of about 21 for cyclic liquefaction occurrence. Table F.3 
presents estimates of equivalent Vs; for blow count of 21. The correlation by Ohta and Goto 
(1978) suggested equivalent Vs; values of 195 m/s for Holocene sands. The stress-corrected 
crosshole compiled by Sykora (1987b) for Holocene sands and non-plastic silty sands below the 
ground water table, with (N,)so between 16 and 26, exhibited an average value of 199 m/s and 
standard deviation of 36 m/s. From these estimates, a Vs; value of 195 m/s is assumed 
equivalent to an (1)60 value of 21 in non-plastic soils with FC > 35 %. 


To permit the CRR-Vs; curves for magnitude 7.5 earthquakes shown in Fig. F.1 to have 
Vs; values between 195 m/s and 210 m/s at CRR near 0.6, values of Y's; are assumed to range 
linearly from 200 m/s to 215 m/s, AS hia f The relationship between )’s, and fines content, 
FC, can be expressed by: 


Vs; = 215 m/s for sands with FC <5 % (F.4a) 
Vy =215 -0.5(FC-5) m/s for sands with 5 % < FC <35% (F.4b) 
Vs, = 200 m/s for sands and silts with FC > 35 % (F.4c) 
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To illustrate how well the recommended CRR-Vs; curves defined by Eqs. (F.1) and (F.4) 
fit the case history data, the data separated by soil type, are presented in Figs. F.4 through F.7. 
The recommended curves provide reasonable bounds for all the case history data above a CSR 
value of 0.35, supporting the use of the suggested 5, values for sands and silts, as well as 


gravels. The use of these )’s, values for gravels is discussed below. 


Table F.3 - Estimates of Equivalent Vs; for Holocene Sands and Gravels Below the Ground 


Water Table with Corrected SPT Blow Count of 21. 


Equivalent Vs; 
Estimate (m/s) 
(3) 


Reference 
(1) 


Relationship 
(2) 


Ohta & Vs = 69 (N°)? 2°)? FLF> 195 
Goto N; = SPT blow count ... for Holocene 
(1978); also measured in Japanese | sands 
given in practice 
report by z =depth, m 214 
Sykora F;, = 1.00 for Holocene- ...for Holocene 
(1987a, age soils gravels 
page 29) F2 = 1.085 for sands; 
1.189 for gravel 
... best-fit relationship for 289 
sets of SPT and Vs 
measurements from Japan 
Sykora Correlation between (N7)60 and 199 
(1987b, crosshole Vs, normalized to .. for Holocene 
page 90); effective overburden stress, sands and non- 
This Report | measurements for Holocene plastic silty sands 
sands and non-plastic silty below the water 
sands below the water table at | table 
sites in U.S.A.; 31 sets of ... Standard 
deviation is 
36 m/s 
Rollinset | Vs= 53 (Noo)”? (6) 217 
al. (1998a) | .. best-fit relationship using ... for Holocene 
equivalent Neo-values from gravels 
Becker Penetration Tests and | --.most of data lie 
Vs measurements; 186 points | within +25 % of 
from 7 Holocene gravel sites relationship 
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Assumptions 
(4) 


1. N;= 60/67 Neo 
2. Neo =21 


Sy 


to 


WN = 


z = 10 m is equivalent 


an overburden stress of 
100 kPa 


. All measurements are 


from below the ground 
water table 


. Average for Vs; values 


with (N))60 between 
16 and 26 


.o', = 100 kPa 


a Neo =21 
.o', = 100 kPa 
. All measurements are 


from below the ground 
water table 
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Fig. F.4 - Curve Recommended for Calculation of CRR from Vs; Measurements in Sands 
and Silts with FC > 35 % along with Case History Data. 
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Data Based on: bd. 
CSR adjusted by dividing by My = 7.5 
MSF = (M,,/7.5)2-96 pape 
My = 5.9 to 8.3 Content (%) 
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sandy soils 
Average values of 
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Fig. F.5 - Curve Recommended for Calculation of CRR from Vs; Measurements in Sands 
with FC = 20 % along with Case History Data. 
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Data Based on: 
CSR adjusted by dividing by 
MSF = (M,,/7.5y2-56 
Mw = 5.9 to 8.3 
Uncemenied, 
Holocene-age 
Sandy soils 
Average values of 
Vs1and amax 


Liquefaction 


No 
Liquefaction 


Field 
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e Liquefaction 
ONo liquefaction 


Content, % 
mo <o 


Cyclic Stress or Resistance Ratio, CSR or CARR 


0 100 200 300 


Overburden Stress-Corrected Shear Wave 
Velocity, Vs1, m/s 


Fig. F.6 - Curve Recommended for Calculation of CRR from Vs; Measurements in Sands 
with FC < 5 % along with Case History Data. 
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Fig. F.7 - Curves Recommended for Calculation of CRR from Vs; Measurements in Gravels 
with FC <5 % and FC = 20 % along with Case History Data. 
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F.1.2 Gravelly Soils 


Although the )’s, values given in Eq. (F.4) were determined for sandy soils, the results 
presented in Fig. F.7 indicate that these limits also represent reasonable limits for gravelly soils 
divided into the same categories based on fines content. This might be considered rather 
surprising, based on the penetration-Vs correlations presented in the literature for gravelly soils. 
For instance, as noted in Table F.1, the correlation by Ohta and Goto (1978) suggested a Vs, 
value of 227 m/s for Holocene gravels at an equivalent (N1)60 of 30. Similarly, the correlation by 
Rollins et al. (1998a) provided a best-fit value of 232 m/s for Holocene gravels. On the other 
hand, all the liquefaction case history data shown in Figs. F.4 through F.7 exhibit Vs, values of 
about 200 m/s or less, suggesting that 230 m/s may be inappropriately high. 


To investigate further the value of )’§, in gravelly soils, laboratory studies involving Vs 
measurements in gravelly soils were reviewed. Kokusho et al. (1995b) clearly showed that the 
shear wave velocity of gravelly soils varies greatly and is highly dependent on the particle 
gradation. Weston (1996) showed similar results for coarse sands with gravels. In both cases, 
the results show that increasing the uniformity coefficient can significantly increase the shear 
wave velocity in medium dense to dense gravels. On the other hand, very loose gravelly soils, 
even well-graded gravels, can exhibit shear wave velocities similar to those of loose sands 
(Kokusho et al., 1995b). The case history data presented in Fig. F.7 supports the premise that 
gravelly soils that are loose enough to exhibit significant liquefaction effects (boils, ground 
cracks, etc.) have shear wave velocities similar to loose sands. Hence, the authors recommend 
the boundaries developed for sandy soils as preliminary boundaries for gravelly soils. However, 
additional work is clearly needed to understand the relationship between Vs, and liquefaction 
resistance of gravels. 


F.2 CURVE FITTING PARAMETERS a AND b 


The curve fitting parameters a and 5 in Eq. (F.1) can be approximated from the case 
history data assuming the values of 5, given in Eq. (F.4) and a MSF relationship. Three MSF 
relationships representing the range of proposed magnitude scaling factors (see Section 2.3.1) 
are considered below to establish the values of a and b in Eq. (F.1). 


F.2.1 Magnitude Scaling Factors Recommended by 1996 NCEER Workshop 
F.2.1.1 Lower Bound of Recommended Range—Figure F.1 presents the case history 
data for magnitude 5.9 to 8.3 earthquakes adjusted using the lower bound for the range of 


magnitude scaling factors recommended by the 1996 NCEER workshop (Youd et al., 1997). 
The lower bound is defined by Eq. (2.9) with n = -2.56, as discussed in Section 2.3.1. Also 
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shown in Fig. F.1 are three recommended CRR-Vs, curves for earthquakes with magnitude near 
7.5 and various fines content. The three curves were determined through an iterative process of 
varying the values of a and 6 until nearly all the liquefaction case histories were bound by the 
curves with the least amount of non-liquefaction case histories in the liquefaction region. The 
final values of a and b used to draw the curves were 0.022 and 2.8, respectively. 


Of the 99 liquefaction case histories, only two liquefaction case histories incorrectly lie in 
the no liquefaction region. The two liquefaction case histories shown in Fig. F.1 that incorrectly 
lie in the no liquefaction region are two sites at Treasure Island (UM0S and UMO9) where 
liquefaction was marginal during the 1989 Loma Prieta earthquake (4, = 7). The sites are 
located along the perimeter of Treasure Island. Mapped liquefaction effects generated by the 
1989 earthquake near the UMOS site are ground cracks with 50 to 90 mm of horizontal 
displacement (R. D. Hryciw, personal communication to R. D. Andrus, 1998; Power et al., 
1998). The nearest mapped sand boil is located 60 m away from the site. At the UMO9 site, as 
much as 90 mm of vertical displacement was observed adjacent to a building located 60 m inland 
from the site. These displacements are small compared to the meters of displacement that are 
expected to occur during larger ground shaking. Thus, liquefaction was marginal at the UMOS 
and UMO9 sites, and sloping ground may have been a factor. It is interesting to note that similar 
incorrect evaluations also are obtained when one uses the SPT and CPT data for these two sites. 
The SPT- and CPT-based evaluations for the UM05 site are discussed in Section F.3.1. 


F.2.1.2 Upper Bound of Recommended Range—Figure F.8 presents the case history 
data for magnitude 5.9 to 8.3 earthquakes adjusted using the upper bound for the range of 
magnitude scaling factors recommended by the 1996 NCEER Workshop (Youd et al., 1997). 
The upper bound is defined by Eq. (2.9) with m = -3.3, as discussed in Section 2.3.1. Also 
shown in Fig. F.8 are the three CRR-Vs, curves from Fig. F.1. Many case histories plot lower in 
Fig. F.8 than in Fig. F.1, because the M,, is less than 7.5 for most of the data. The downward 
shift in the liquefaction data points near the curves at CRR of about 0.08 is less than 0.01. This 
difference is not significant, and is within the accuracy of the plotted case history data. 


F.2.2 Revised Magnitude Scaling Factors Proposed by Idriss (1999) 


Figure F.9 presents the case history data for magnitude 5.9 to 8.3 earthquakes adjusted 
using the revised magnitude scaling factors and stress reduction coefficients proposed by Idriss 
(1999, as discussed in Section 2.3.1. Also shown in Fig. F.9 are the three CRR-Vs, curves from 
Fig. F.1. Many of the case history data shown in Fig. F.9 plot higher than case history data in 
Fig. F.1, because the M, is less than 7.5 for most of the data) The upward shift in the 
liquefaction data points near the curves at CRR of about 0.08 is less than 0.01. Again, this 
difference is not significant, and is within the accuracy of the plotted case history data. 
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Fig. F.8 - Curves Recommended for Calculation of CRR from Shear Wave Velocity 
Measurements Along with Case History Data Based on Upper Bound Values 
of MSF for the Range Recommended by the 1996 NCEER Workshop (Youd 
et al., 1997) and rg Developed by Seed and Idriss (1971). 
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Fig. F.9 - Curves Recommended for Calculation of CRR from Shear Wave Velocity 
Measurements Along with Case History Data Based on Revised Values of MSF 
and rg Proposed by Idriss (1999). 
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F.2.3 Comparison of Magnitude Scaling Factors 


The proposed relationships for MSF can be compared directly by combining them with 
the appropriate stress reduction coefficient into one factor. This factor is the product of rz and 
the reciprocal of MSF. Figure F.10 presents values of ra/MSF for the range recommended by 
the 1996 NCEER Workshop (Youd et al., 1997) and those proposed by Idriss (1999). As 
shown in the figure, there is not much difference between the two sets of r4/MSF values for 
magnitude of 7.5 and depth less than 11 m. At magnitudes near 5.5 and shallow depths, the 
difference between r./MSF values proposed by Idriss (1999) and values recommended by the 
1996 NCEER Workshop is as much as 50 %. Thus, at magnitudes less than about 7, the 
difference in using values of MSF and rz proposed by Idriss (1999) and those adopted by the 
NCEER Workshop (Youd et al., 1997) is significant in the calculation of CSR. 


For example, Fig. F.11 presents two liquefaction resistance curves for earthquakes with 
magnitude near 5.5 and clean soils (FC < 5 %). The upper curve was obtained by multiplying 
values of CRR defining the curve for FC < 5% in Fig. F.1 by 2.2, the lower MSF recommended 
by the 1996 NCEER Workshop for magnitude 5.5 earthquakes (see Eq. (2.9) with n = -2.56). 
The lower curve was obtained by multiplying values of CRR defining the curve for FC < 5 % in 
Fig. F.1 by 1.68, the MSF proposed by Idriss (1999) for magnitude 5.5 earthquakes (see Eq. 
(2.10)). Also shown in Fig. F.11 are the available case history data for clean sands determined 
using average stress reduction coefficients proposed by Seed and Idriss (1971) and Idriss 
(1998). The two curves in Fig. F.11 exhibit differences in CRR of about 0.02 at Vs; = 100 m/s 
and 0.1 at Vs; = 200 m/s. 


F.3 RECOMMENDED CRR-Vs; CURVES 


From the discussion presented above, the recommended CRR-Vs; curves are defined by 
Eqs. (F.1), (2.9) and (F.4) with a = 0.022, b = 2.8, and n= -2.56. The recommended curves for 
moment magnitudes ranging from 5.5 though 8 are presented in Figs. F.12 through F.17, 
respectively, along with the case history data. The value of -2.56 for m is recommended for 
determining magnitude scaling factors because it provides more conservative CRR-Vs; curves 
than -3.3, which is the 7 value defining the upper bound of the range of MSFs suggested by the 
1996 NCEER Workshop (Youd et al., 1997) for magnitudes less than 7.5. Although the 
magnitude scaling factors defined by Eq. (2.9) with n = -2.56 provide less conservative CRR-Vs; 
curves than the factors proposed by Idriss (1999) for magnitudes less than 7.5, the factors 
determined of Ambraseys (1988), I. M. Idriss (personal communication to T. L. Youd, 1995), 
Arango (1996), Youd and Noble (1997), and Andrus and Stokoe (1997; as indicated by the very 
conservative CRR-Vs; curves shown in Figs. F.13 and F.14) supported their use. 
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Fig. F.10 - Variation of r,/MSF with Depth for Various Magnitudes and Proposed Relationships. 
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Fig. F.11 - Comparison of Liquefaction Resistance Curves and Case History Data 
for Procedures Recommended by the 1996 NCEER Workshop 
(Youd et al., 1997) and the Revised Procedures Proposed by Idriss (1999) 
for Clean Sands and Earthquakes with Magnitude Near 5.5. 
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Fig. F.12 - Case History Data for Earthquakes with Magnitude Near 5.5 Based on 
Overburden Stress-Corrected Shear Wave Velocity and Cyclic Stress Ratio 
with Recommended Liquefaction Resistance Curves. 
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Fig. F.13 - Case History Data for Earthquakes with Magnitude Near 6 Based on 
Overburden Stress-Corrected Shear Wave Velocity and Cyclic Stress Ratio 
with Recommended Liquefaction Resistance Curves. 
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Fig. F.14 - Case History Data for Earthquakes with Magnitude Near 6.5 Based on 
Overburden Stress-Corrected Shear Wave Velocity and Cyclic Stress Ratio 
with Recommended Liquefaction Resistance Curves. 
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Fig. F.15 - Case History Data For Earthquakes with Magnitude Near 7 Based on 
Overburden Stress-Corrected Shear Wave Velocity and Cyclic Stress Ratio 
with Recommended Liquefaction Resistance Curves. 
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Fig. F.16 - Case History Data for Earthquakes with Magnitude Near 7.5 Based on 
Overburden Stress-Corrected Shear Wave Velocity and Cyclic Stress Ratio 
with Recommended Liquefaction Resistance Curves. 
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Fig. F.17 - Case History Data for Earthquakes with Magnitude Near 8 Based on 
Overburden Stress-Corrected Shear Wave Velocity and Cyclic Stress Ratio 
with Recommended Liquefaction Resistance Curves. 
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The recommended CRR-Vs; curves shown in Figs. F.12 through F.17 are dashed above 
CRR of about 0.35 to indicate that they are based on limited field performance data. The curves 
do not extend much below 100 m/s, because there are no field data to support extending them to 
the origin. It is important to note that these boundary curves are for extreme behavior, where 
boils and ground cracks occur. 


F.4 CORRELATIONS BETWEEN Vs; AND PENETRATION RESISTANCE 


One can obtain correlations between Vs, and corrected penetration resistance from the 
recommended CRR-Vs; relationships given in Fig. F.1 and 1996 NCEER Workshop (Youd et 
al., 1997) recommended SPT- and CPT-based relationships for magnitude 7.5 earthquakes by 
plotting values with equal CRR. 


F.4.1 Corrected SPT Blow Count 


Figure F.18 presents the correlation of Vs; with (N1)60 for clean soils (< 5 % fines), based 
on the recommended CRR-Vs; and CRR-(N,)60 relationships. Also shown are the field data and 
mean curve for sands with less than 10 % non-plastic fines from Fig. F.2. The correlation 
derived from the CRR relationships lies between the mean and the mean +15,.; curves. The 
flatter slope below (Nj)6o of 6 exhibited by the CRR-based correlation can be explained by 
different assumed minimal values of CRR. The CRR-Vs; relationship for magnitude 7.5 
earthquakes and FC < 5 % shown in Fig. F.1 provides a CRR of 0.033 for Vs; = 100 m/s, the 
lowest Vs; value shown in the figure. The 1996 NCEER Workshop recommended a CRR value 
of 0.05 for (Vi)6o = 0. The difference between minimal values of CRR is small, and is near the 
accuracy of both procedures. 


The CRR-based correlation shown in Fig. F.18, along with the plotted field data, provide 
a simple method of comparing the Vs;- and (N,)6o-based liquefaction evaluation procedures. 
Both procedures provide similar predictions of liquefaction potential, when the data point lies on 
the CRR-based curve. When the data point plots below the CRR-based curve, the Vs;-based 
liquefaction evaluation procedure provides the more conservative prediction. When the data 
point plots above the CRR-based curve, the SPT-based liquefaction evaluation procedure 
provides the more conservative prediction. Because most of the data points shown in Fig. F.18 
plot below the CRR-based curve, the Vs;-based procedure provides an overall more conservative 
prediction of liquefaction resistance than does the SPT-based procedure for these sites. 
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Fig. F.18 - Relationships Between (N;)59 and V;; for Clean Sands Implied by the Recommended 
CRR-Vg, Relationship and the 1996 NCEER Workshop Recommended CRR-(N,)¢0 


Relationship (Youd et al., 1997) with Field Data for Sands with Less than 10 % 
Fines. 
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The data point for the Treasure Island UM0S site, which incorrectly lies in the region of 
no liquefaction shown in Fig. F.1, plots just below the CRR-based curve, as shown in Fig. F.18. 
Thus, this case history also incorrectly plots in the region of no liquefaction on the SPT-based 
liquefaction evaluation chart. Furthermore, the SPT-based procedure provides a slightly less 
conservative prediction of liquefaction resistance than the shear-wave-based procedure for this 
case history. 


Although the CRR-based curve shown in Fig. F.18 generally trends parallel to the mean 
curve, there is a small hump between corrected blow counts of 8 and 26. This hump suggests - 
that either the CRR-Vs, relationship is more conservative or the CRR-(N1)¢o0 relationship is less 
conservative in this range. 


Similarly, as shown in Fig. F.19, a correlation between Vs; with (V1)60 for soils with > 35 
% fines can be derived from the recommended CRR-Vs; and CRR-(N;)oo relationships. Figure 
F.19 provides the basis for the method of estimating the cementation and aging correction 
factor, K., suggested in Section 2.3.4 (see Fig. 2.6). 


F.4.2 Normalized Cone Tip Resistance 


Figure F.20 presents the correlation of Vs; with g.:~ for clean sands with median grain 
_ size, Dso, between 0.25 mm and 2.0 mm, based on the recommended CRR-Vs,; and CRR-q-in 
relationships. Also shown are the field data and mean curve for clean sands with less than 10 % 
non-plastic fines from Fig. F.3. The correlation derived from the CRR relationships lies between 
the mean and the mean +15,., curves for Vs; = 170 m/s, indicating that the Vs;-based procedure 
provides an overall more conservative prediction of liquefaction resistance than does the CPT- 
based procedure for these sites. For Vs; < 170, the CRR-based correlation lies close to the mean 
curve, indicating that both procedure provide an overall similar prediction. The slope of the 
CRR-based correlation below q.iv of 20 may be explained by the different assumed minimal 
values of CRR, as discussed in Section F.4.1. 


The data point for the Treasure Island UM05 site, which incorrectly lies in the region of 
no liquefaction shown in Fig. F.1, plots on the CRR-based curve, as shown in Fig. F.20. Thus, 
the CPT-based procedure provides a similar incorrect prediction of no liquefaction for this site. 


Similarly, as shown in Fig. F.21, a correlation between Vs; with q.1w for soils with > 35 
% fines can be derived from the recommended CRR-Vs; and CRR-(N,)60 relationships. Figure 
F.21 provides the basis for the method of estimating the cementation and aging correction 
factor, K., suggested in Section 2.3.4 (see Fig. 2.7). 
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Fig. F.19 - Relationships Between (N;)<9 and V;, Implied by the Recommended CRR-V., 
Relationship and the 1996 NCEER Workshop Recommended CRR-(N})¢60 
Relationship (Youd et al., 1997). 
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Fig. F.20 - Relationships Between q,), and V;, for Clean Sands Implied by the Recommended 
CRR-Vs; Relationship and the 1996 NCEER Workshop Recommended CRR-q,.y 
Relationship (Youd et al., 1997) with Field Data for Sands with Less than 10 % 
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Fig. F.21 - Relationships Between q,,, and V;; Implied by the Recommended CRR-V, 
Relationship and the 1996 NCEER Workshop Recommended CRR-q,;y 
Relationship (Youd et al., 1997). 
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F.5 SUMMARY 


The development of the recommended CRR-Vs; curves was outlined in this appendix. 
The recommended curves are based on a modified relationship between shear wave velocity and 
cyclic stress ratio for constant average cyclic shear strain suggested by R. Dobry. They are 
defined by Eqs. (F.1), (F.4) and (2.9) with a = 0.022, b = 2.8, and n = -2.56. The curve fitting 
parameters a and Bb are determined through an iterative process that involved varying their 
values until nearly all the liquefaction case histories were bound by the curves with the least 
amount of non-liquefaction case histories in the liquefaction region. Three MSF relationships 
are considered in determining the values of a and b. Equation (F.4), which provides a 
relationship between the limiting upper Vs; value and fines content, is based, in part, on the case 
history data and, in part, on penetration shear wave velocity correlations. From penetration-Vs 
correlations, the recommended CRR-Vs, curves appear to be somewhat more conservative than 
the penetration-based curves recommended by the 1996 NCEER Workshop (Youd et al., 1997). 
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APPENDIX G 


PROBABILITY-BASED LIQUEFACTION EVALUATION 


This appendix presents three probability models for the case history data listed in 
Appendix H. The probability models are based on the work of Juang et al. (2001a; 2002), and 
are derived using logistic regression and Bayesian interpretation techniques. They are compared 
with the deterministic evaluation curve by Andnmus et al. (1999) for clean soils (FC < 5 %) shown 
in Fig. 2.3. The probability models provide a means of objectively calibrating the deterministic 
liquefaction evaluation curve. 


To develop the probability models, values of Vs; are adjusted to a clean soil equivalent. 
The procedure for adjusting Vs; values involves two steps. First, a CRR value is determined 
using Eq. (2.8) for each case history. Second, for each value of CRR, a clean soil equivalent Vs; 
value is determined using Eq. (2.8) with /’;, = 215 m/s. Thus, this adjustment procedure, 
maintains the ratio of CRR to CSR (or factor of safety). The adjustment procedure can be 
expressed by: 


Vsics = Ky Vs: (G.1) 
where 


Vsics = the equivalent clean soil value of Vs;, and 
Kj; | =a fines content correction to adjust Vs; values to a clean soil equivalent. 


Values of Ky. can be approximated using the following equation (Juang et al., 2001a; 2002): 


Ky ="1, for FCS % (G.2a) 
Ky = 1+ (FC - 5) fVs,), for FC = 6 % to 34 % (G.2b) 
Ky = 1+ 30 fKVs,), for FC > 35% (G.2c) 
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where 
V. eve 
= 0.009 - 0.0109 | -S4+ |+ 0.0038 | —S/ G3 
KY si) [F. [Fa Soe 


Equations (G.1) through (G.3) provide an approximate mathematical description of the 
adjustment procedure. The adjusted case history data are plotted in Figs. G.1 and G.2 along 
with two probability models determined using logistic regression. The logistic regression-based 
probability models, as well as a Bayesian-based probability model, are discussed below. 


G.1 LOGISTIC REGRESSION MODELS 
G.1.1 Logistic Regression Model 1 
The first logistic regression-based probability model, called Model 1, is similar in form to 


the model used by Liao et al. (1988) for analyzing SPT-based case histories. The probability 
equation for Model 1 is given by (Juang et al., 2001a; 2002): © 


In Fi =@Q) ot Qa? Vs1.cs ts Q3 In(CSR7.5) (G.4) 
1 cp? P, 
where 
P, = the probability that liquefaction will occur, 
Qa}, Az, a3 = regression coefficients, and 
CSR7.5 = CSR adjusted to M,, = 7.5. 


The mean values of a;, a2, and a3 are 14.8967, -0.0611, and 2.6418, respectively. The standard 
deviations associated with the coefficients are 2.1637, 0.0098, and 0.4268, respectively. The 
Nagelkerke coefficient (equivalent to R*) of this regression is 0.58. Probability curves for 
Model 1 are presented in Fig. G.1. From the figure, Model 1 appears to provide reasonable P;, 
curves within the limits of most of the data. However, the Model 1 curves may be 
inappropriately too conservative at high values of Vs;,cs (say > 200 m/s), since a corrected 
velocity of 210 m/s is considered equivalent to a corrected blow count of 30 in clean sands and 
liquefaction is generally assumed not possible above this value. 
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Fig. G.1 - Logistic Regression Model 1 and Case History Data Adjusted for Fines Content. 
(after Juang et al., 2001a) 


G.1.2 Logistic Regression Model 2 


To investigate the influence that the form of a regression equation might have on P, 
curves, the analysis is repeated using a slightly different equation. The probability equation for 
this second logistic regression model, called Model 2, is defined by (as suggested by William 
Guthrie, NIST, to R. D. Andrus, June 1998): 


| s iF by + by Vsics + bs In(CSR75) + by [In(CSR75)]° (G.5) 


Leet 8G 
where 
b;, b2, b3, b4 += regression coefficients. 


The mean values of 5;, b2, b; and b, determined by Juang et al. (2001a; 2002) are 10.0155, - 
0.0643, -3.9534, and -1.8381, respectively. The standard deviations associated with the 
coefficients are 2.6102, 0.0107, 2.1738, and 0.6302, respectively. The Nagelkerke coefficient of 
this regression is 0.61. A coefficient of 0.61 is slightly greater than 0.58, suggesting a slightly 
stronger correlation for Model 2 than Model 1. Figure G.2 presents P, curves defined by Eq. 
(G.5). The Model 2 curves exhibit steeper-slopes than Model 1 curves above a CSR value of 
about 0.1. They reach a maximum Vs;,cs value at CSR of about 0.33. Above CSR of 0.33, the 
curves trend to the left, decreasing in Vs;cs with increasing CSR. Nevertheless, the results 
clearly show that P, curves determined by logistic regression depend on the form of the 
regression equation. While Model 2 provides another possible probability model, one would 
expect P, curves to slope towards higher values of Vs;,cs with increasing CSR rather than extend 
vertically, as suggested by the dashed lines in Fig. G.2. 


G.2 BAYESIAN MAPPING MODEL 


Juang et al. (1999) pioneered a Bayesian interpretation approach for mapping factor of 
safety, Fs, to Pr. In their approach, values of Fs are first determined for the liquefaction and 
non-liquefaction case histories using a deterministic evaluation curve. The Vs-based curve 
shown in Fig. 2.3 is the deterministic curve used in this case. Values of P, are then estimated 
from the probability density functions of F's for liquefaction and non-liquefaction case histories 
using Bayes’ theorem. With the assumption of equal prior probability, the P,-Fs mapping 
function can be expressed as (Chen and Juang, 2000; Juang et al., 2000a): 


epee Dist 


Tepes Tee G.6 
t,(Fs)+ fa Fs) E59) 


Liquefaction 


Cyclic Stress or Resistance Ratio, CSR or CRRA 


100 150 200 250 300 
Modified Shear Wave Velocity, Vs1,cs, m/s 


Fig. G.2 - Logistic Regression Model 2 and Case History Data Adjusted for Fines Content. 
(after Juang et al., 2001a) 
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where 


fi(Fs) = the probability density function of the calculated F's for the liquefaction case 
histories, and 

Jut(F's) = the probability density function of the calculated F's for the non-liquefaction 
case histories. 


An analysis of the 225 case histories yields the probability density functions shown in 
Figs. G.3a and G.3b for the liquefaction and non-liquefaction cases, respectively. Applying Eq. 
(G.6), the predicted probability of liquefaction for each case history is obtained. Values of P, 
and Fs for the case histories are plotted in Fig. G4. The relationship formed by the Fs -P, 
values can be approximated by (modified from Juang et al., 2001a): 


[9h oe (G7) 


In Eq. (G.7), a F's value of 1 corresponds to points on the deterministic curve. Thus, on 
average, the Andrus et al. (1999) curve for clean soils (see Fig. 2.3) 1s characterized with a P; 
value of 26 % based on Eq. (G.7). The value of 26 % is slightly less than 30 % initially 
determined by Juang et al. (2001a). Subsequent analysis (Juang et al., 2002) revealed that a few 
of the calculated P; values corresponding to low Fs values (see Fig. G.4) were unreasonably 
influencing the coefficients given in Eq. G.7. Thus, Eq. G.7 has been modified slightly from the 
preliminary equation proposed by Juang et al. (2001a). 


Figure G.5 compares the Fs-P; relationship defined by Eq. (G.7) for the Vs-based 
recommended curve (Andrus et al., 1999) with the F's-P,; relationship developed by Juang et al. 
(2000) for the SPT-based recommended curve (Seed et al., 1985; Youd et al., 2001). There is 
remarkable agreement between the Vs- and SPT-based relationships. From Fig. G.5, the SPT- 
based recommended curve is characterized with an average P; value of 31 %. These findings 
suggest that the Vs-based deterministic evaluation curves are somewhat more conservative than 
the SPT-based curves. 


Equation (G.7) provides an important link between the probabilistic and deterministic 
methods. One can obtain a family of P; curves for probability-based design by combining Eqs. 
(2.8), (2.12) and (G.7). The family of P, curves for magnitude 7.5 earthquakes and soils with 
FC <5 % is presented in Fig. G.6. These curves, called the Bayesian Mapping Model, slope to 
the right with increasing CSR, which seems reasonable. They converge to a Vs, value of 215 
m/s, the assumed value of Vs, for clean soils, at high values of CSR. 
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Fig. G.3 — Probability Density Functions and Calculated Factor of Safety Distributions for the 
225 Case Histories. 
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Fig. G.4 — Calculated Values of Fs and P, for the 225 Case Histories with Best-Fit Relationship 
(Juang et al., 2002). Note that Fs is based on the Recommended Deterministic 
Evaluation Curves Shown in Fig. 2.3, and P, is based on Eq. (G.6). 
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Fig. G.5 - Relationship Between P, and F's for the Vs-based Procedure (Juang et al., 2002) and 
SPT-Based Procedure (Juang et al., 2000a) Determined Using Bayes’ Theorem. 
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Fig. G.6 - Bayesian Mapping Model and Case History Data Adjusted for Fines Content. (Juang 
et al., 2002) 
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G.3 COMPARISON OF PROBABILITY MODELS 


Figure G.7 compares the logistic regression model curves for P, = 26 % with the 
Bayesian Mapping Model curve for P; = 26 %, which corresponds to the deterministic curve 
developed by Andrus et al. (1999) for soils with FC < 5%. The Bayesian Mapping Model curve 
lies between the two logistic regression curves below a Vs;-; value of about 195 m/s, indicating 
close agreement between the three probability models. Above 195 m/s, the Bayesian Mapping 
Model curve closely follows the logistic regression Model 2 curve. Thus, the logistic regression 
models support the Bayesian Mapping Model in characterizing the deterministic curve proposed 
by Andrus et al. (1999) and Andrus and Stokoe (2000) as a 26 % probability of liquefaction 
curve. 


The tendency for the P; curves to converge to some limiting upper value reflects the 
tendency of dense soils to exhibit dilative behavior at large strains, causing negative pore-water 
pressures. It seems reasonable that the P; curves should not continue to diverge with increasing 
Vs, Or penetration resistance, but should converge somewhat to reflect the behavior of dense 
soils, as suggested by the curves shown in Fig. G.6. The wider spread exhibited in logistic 
regression-based P; curves at high values of Vs and CSR is believe to be the result of an inherent 
property of these models, and not a real-world phenomenon. Thus, the Bayesian Mapping 
Model (Fig. G.6) is considered to be an improvement over the logistic regression models, and is 
suggested for engineering risk-based design. 
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Fig. G.7 - Comparison of the Probability-Based Logistic Regression and Bayesian Mapping 
Models for P; = 26 % with the Deterministic Curve Developed by Andrus et al. 
(1999) for Clean Soils. (modified from Juang et al., 2001a) 
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APPENDIX H 


SUMMARY OF CASE HISTORY DATA 


~ 


Table H.1 presents a summary of case history data described in Appendix E, and used in 
Appendix F to establish the recommended liquefaction resistance curves. This database is 
expanded and modified from the database presented by Andrus and Stokoe (1997). Most of the 
modifications are minor with the intent to have the data conform to the guidelines presented in 
this document. The major modifications are based on new information or correction of an error 
in calculations. Some case histories included in the earlier database by Andrus and Stokoe have 
been omitted due to one of the three following reasons: (1) The reported average downhole Vs 
measurement is for a depth interval much greater than the identified critical layer. (2) The 
critical layer is likely older than 10 000 years and contains carbonate. (3) The location of the 
critical layer or field behavior is uncertain. The case history data presented in Table H.1 are 
essentially the same as the data presented in the draft guidelines (Andrus et al., 1999), with only 
a few minor changes. References for the case history data are given in Table E.1. 
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istories. 


Table H.1 - Summary Information for Vs -Based Liquefaction and Non-Liquefaction Case H 
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